REGULARITY OF SOLUTIONS TO THE POLYHARMONIC EQUATION IN 

GENERAL DOMAINS 



SVITLANA MAYBORODA AND VLADIMIR MAZ'YA 

Abstract. The present paper establishes boundedness of [m - | + 5] derivatives for the solutions 
to the polyharmonic equation of order 2m in arbitrary bounded open sets of R", 2 < n < 2m + 1, 
without any restrictions on the geometry of the underlying domain. It is shown that this result is 
sharp and cannot be improved in general domains. Moreover, it is accompanied by sharp estimates 
on the polyharmonic Green function. 
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1. Introduction 

Higher order differential operators are important in physics and in engineering [34] and have 
been integrated in many areas of mathematics, including conformal geometry (Paneitz operator, 
2-curvature [8], [9]), free boundary problems [1], non-linear elasticity [50], [10], [5], to mention 
just a few. Unfortunately, in spite of evident demand, the properties of higher order PDEs on gen- 
eral domains remained largely beyond reach. Their investigation brought challenging hypotheses 
and surprising counterexamples, and few general positive results. For instance, Hadamard's 1908 
conjecture regarding positivity of the biharmonic Green function [16] was actually refuted in 1949 
(see [14], [15], [45]), and later on the weak maximum principle was proved to fail as well, at least 
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in high dimensions [32], [44]. Another curious feature is a paradox of passage to the limit for 
solutions under approximation of a smooth domain by polygons [6], [28]. 

For a long time, almost all results for higher order PDEs pertained to piecewise smooth domains 
(see, e.g., [19], [20], [33]). The recent two decades have witnessed a great burst of activity in the 
study of boundary value problems on Lipschitz domains, which brought some powerful and deli- 
cate estimates for the bi-Laplacian and for the polyharmonic operator due to B. Dahlberg, C. Kenig, 
J.Pipher, G. Verchota, Z. Shen, M.Mitrea, and others ([11], [53], [40], [41], [42], [17], [44], [43], 
[54], [12], [2], [55], [46], [47], [48], [49], [35], [27], [56], [18], [36]). Unfortunately, none of these 
results could be extended to domains of more complicated geometry. 

The present paper establishes sharp pointwise estimates on the solutions to the polyharmonic 
equation and their derivatives in arbitrary bounded open sets, without any restrictions on the 
geometry of the underlying domain. It is shown that these estimates are sharp and can not be 
improved. 

Let us recall the second order case. One of the fundamental results of elliptic theory is the 
maximum principle for harmonic functions. It holds in arbitrary domains and guarantees that 
every solution to the Dirichlet problem for the Laplace's equation, with bounded data, is bounded. 
A similar statement is valid for a broad class of second order elliptic operators. 

In 1960 an appropriate, weak, version of the maximum principle was obtained for higher order 
equations on smooth domains ([3], see also [37], [38]). Roughly speaking, it amounts to the 
estimate 

(1.1) \\V"- i u\\ L ^ ) <C\\V n - 1 u\\ L ^ m) , 

where u is a solution of an elliptic differential equation of order 2m such that |V" ,_1 i/| is continuous 
up to the boundary, V m_1 w = {d a u}\ a \ =m -\ is a vector of all partial derivatives of u of order m - 1 and 
we adopt the usual convention that the zeroth order derivative of u is u itself. The formulation in 
non-smooth domains is somewhat trickier (see, e.g., [39]). However, in any setting, the weak max- 
imum principle would always guarantee that the solution with "nice" data has bounded derivatives 
of order m - 1 . In the early 1 990s, (1.1) was extended to three-dimensional domains diffeomorphic 
to a polyhedron ([20], [31]) or having a Lipschitz boundary ([42], [44]). In general domains, no 
direct analog of the maximum principle exists (see Problem 4.3, p. 275, in Necas's book [39]). 

Moreover, it turns out that for every elliptic operator of order greater than two the maximum 
principle can be violated, in a four-dimensional cone ([32], see also [29], [41], [44]). It has been 
shown, in particular, that in dimensions n > 4 there are solutions to the polyharmonic equation with 
unbounded (m - l)-st derivatives (cf. (1.1)). This phenomenon raises two fundamental questions: 
whether the boundedness of the (m - l)-st derivatives remains valid in dimensions n < 3, and 
whether there are some other, possibly lower-order, estimates that characterize the solutions when 
n > 4. The main result of the present paper is as follows. 

Theorem 1.1. Let Vibe a bounded domain in R", 2 < n < 2m + 1, and 



(1.2) (— A) m w = / in Q, fe C"(fi), u e W m2 {£l). 
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Then the solution to the boundary value problem (1.2) satisfies 

(1.3) V'"-" /2+1/2 M 6 L°°(Q) when n is odd and V m ' n/2 u e L°°(Q) when n is even. 
In particular, 

(1.4) V m-1 K e L°°(Q) whenn = 2,3. 

Here the space W m ' 2 (Cl), is, as usually, a completion of C^(Q.) in the norm given by ||w||^m.2 (n) = 
||V m w|| L 2 (n) . We note that W" l ' 2 (Q.) embeds into C k (Q.) only when k is strictly smaller than m - f , 
n < 2m. Thus, whether the dimension is even or odd, Theorem 1.1 gains one derivative over the 
outcome of Sobolev embedding. 

The results of Theorem 1 . 1 are sharp, in the sense that the solutions do not exhibit higher smooth- 
ness than warranted by (1.3)— (1.4) in general domains. Indeed, assume that n e [3, 2m + 1] n N is 
odd and let Q c W be the punctured unit ball B t \ {O}, where B r = {x e W : \x\ < r). Consider a 
function 77 6 C ( ^°(5 1 / 2 ) such that 77 = 1 on fi 1/4 . Then let 

(1.5) u(x) := j 1 {x)dT' i ~H\x\ 2m - n ), x e B 1 \ {O}, 

where d x stands for a derivative in the direction of x, for some i = l,...,n. It is straightforward to 
check that u e W" h2 (Q.) and (—A) m u e C~(Q). While l m -'i^-u is bounded, the derivatives of the 
order m - | + | are not, and moreover, V m ~3+i M is not continuous at the origin. Therefore, the 
estimates (1 .3)— (1 .4) are optimal in general domains. 

As for the case when n is even, the results in [20, Section 10.4] demonstrate that in an exterior 
of a ray there is an m-harmonic function behaving as \x\ n ''^^. Thus, upon truncation by the 
aforementioned cut-off 77, one obtains a solution to (1.2) in B\ \ {x\ = 0, ...,x„_i = 0, < x n < 1}, 
whose derivatives of order m - \ + 1 are not bounded. Therefore, in even dimensions (1.3) is a 
sharp property as well. 

Furthermore, Theorem 1.1 has several quantitative manifestations, providing specific estimates 
on the solutions to (1.2). Most importantly, we establish sharp pointwise estimates on Green's 
function of the polyharmonic operator and its derivatives, once again without any restrictions on 
the geometry of the domain. The full list of the estimates is quite extended. For the purposes of 
the introduction, let us just highlight the highest order case. 

As customary, we denote by G(x,y), x, y 6 Q, Green's function for the polyharmonic equation 
and by T its fundamental solution, so that, in particular, G(x, y) - T(x - y), x,y e £1, is the regular 
part of the Green function. By definition, for every fixed y 6 Q the function G(-,y) satisfies 

(1.6) (-A x ) m G(x,y) = 5(x -y), xe Q, 

in the space W m ' 2 (Q,). Here A x stands for the Laplacian in x variable, and similarly we use the 
notation A y , V v , for the Laplacian and gradient in y, and gradient in x, respectively. By d(x) we 
denote the distance from x e Q. to dQ,. 

Theorem 1.2. Let Q c R" be an arbitrary bounded domain. Ifn 6 [3, 2m + 1] n N is odd then 

(1.7) Vr H V;- f+ W,y)-rC*-y)) < — ( : % : , for every x, y e Q, 

max{a(;c), d(y),\x-y\} 
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and, in particular, 



(1.8) 




V 



c 



for every x,y e Q.. 



G(x,y) < 



\x-y\' 



Ifn 



€ [2, 2m] n N is even, then 



(1.9) 



Vr 5 V;- 5 (G(*,y) - Y(x - y)) < C log 1 + 



max{rf(x), d(y), |jc - y|} 



diamQ 



) 



, for every x, y € Q, 



and 



(1.10) 




G(x,y) <Clog 1 + 



min{ J(x), d(y)} 
\x-y\ 



) 



for every x, y e Q. 



The constant C in (1 .7)— (1 . 10) depends on m and n only. In particular, it does not depend on the 
size or the geometry of the domain Q. 

We mention that the pointwise bounds on the absolute value of Green's function itself have been 
treated previously in dimensions 2m + 1 and 2m + 2 for m > 2 and dimensions 5, 6, 7 for m = 2 in 
[25, Section 10] (see also [24]). Also, as for the solutions, there exist results in smooth domains 
[13], [21], [51], [52], in conical domains [30], [20], and in polyhedra [31]. The estimates on the 
derivatives of Green's function in arbitrary bounded domains, provided by Theorem 1.2, are new. 

Furthermore, using standard techniques, the Green's function estimates can be employed to 
establish the bounds on the solution to (1.2) for general classes of data /, such as U for a certain 
range of p, Lorentz spaces etc. We defer the detailed discussion of those to the body of the paper. 

Here we just would like to point out that until recently, the properties stated in Theorems 1.1- 
1.2 seemed accessible exclusively under heavy restrictions on the geometry of the domain (see the 
references above). The present paper rests on a new method, based on some intricate weighted 
integral identities. The biggest challenge, and the core idea, is the proper choice of the weight 
function w. It is very subtle and finely tuned to the underlying elliptic operator in such a way 
that the positivity, or, rather, suitable bounds from below, could be obtained for expressions akin 
to J a (-A) m u(x)u(x) w(x) dx. When such integrals have been considered before (in particular, we 
have established the three-dimensional biharmonic version of the results of this paper earlier in 
[22]), the difficulties of handling general m and n seemed insurmountable. One of the main tech- 
nical achievement of this paper is the novel systematic construction of the weight leading to sharp 
bounds for the solutions. It invokes numerous new aspects: employing induction in eigenvalues 
of the Laplace-Beltrami operator on the sphere, preservation of some positivity properties under 
a change of underlying higher order operator, exploiting delicate peculiarities of (-A) m depend- 
ing on the parity of m,n,m - n/2, and others. The construction appears for the first time in this 
manuscript and is likely to be applicable to general classes of elliptic equations. 

Acknowledgements. We are greatly indebted to Marcel Filoche for the idea relating certain 
positivity properties of one-dimensional differential operators to particular configurations of the 
roots of associated polynomials. It has been reflected in Section 3 and it has ultimately significantly 
influenced our technique. 
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2. Integral inequalities and global estimate: the case of odd dimension. Part I: power weight 

Let us start with a list of notation and conventions used throughout the paper. 

For any domain OcR"a function u e C~(Q) can be extended by zero to W and we will write 

o 

u € CT(K") whenever convenient. Similarly, the functions in W£(Q), m e N, will be extended by 
zero and treated as functions on W or other open sets containing Q without further comments. 

The symbols B r (Q) and S r (Q) denote, respectively, the ball and the sphere with radius r centered 
at Q and C nR (Q) = B R (Q) \ B r (Q). When the center is at the origin, we write B r in place of B r {0), 
and similarly S r := S r (0) and C r ^ R := C ryR (0). 

Let (r, to) be spherical coordinates in R", n > 2, i.e. r = \x\ e (0, oo) and to = x/\x\ is a point 
of the unit sphere In fact, it will be more convenient to use e~', t e R, in place of r, so that 
t = log r" 1 = log W^ 1 . Then by x we denote the mapping 

(2.1) R" 3 x -^U (t,io) e RxS n ~\ n>2. 

The symbols 8 U and V w refer, respectively, to the Laplace-Beltrami operator and the gradient on 

Finally, by C, c, C, and c,-, i e N, we generally denote some constants, possibly depending on 
the order of operator m and the dimension n but not on any other variables and not on the domain, 
unless explicitly stated otherwise. Their exact values are of no importance and can change from 
line to line. Also, we write A « B, if C 1 A < B < C A for some C > 0. 

Theorem 2.1. Assume that m e N and n 6 [3, 2m + 1] n N is odd. Let Qbe a bounded domain in 
R n , O 6 W \ Q, u 6 anJ v = <?( m -H>( M o x^ 1 ). Then 

f (-A) ra M (j) hWW'Ji 

(dfv) dcodt + C v |"| (-6 W - p(p + n-2))vda>dt, 

>.----• nl JrJS"~ 1 n _„n,3 

p- 2 f 2 

where C > is some constant depending on m and n only. 

Proof. Step I. In the system of coordinates (t, to) the polyharmonic operator can be written as 

m-1 

(2.3) (-A) m = (-1) V mf - 2j)(-d t -2j + n-2) + dj). 
Then 

(2.4) I (-A) m M (x) wWlxl -1 = I I £ m ' n (d t ,6 0J )v(t,a))v(t,a))da)dt, 
Jr" Jr Js"" 1 

with 

m-1 / -. ^ v 

(2.5) £ m '"(«9,, dj = (-l) m f] i(-d t + m -l + -- 2j)(-d t + m + l---2j) + sA 

y=0 ' ' 
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Denote by v p ; the coefficients of the expansion of v into spherical harmonics: 

oo p 

(2.6) v(t, co) = £ J] V(0if t 6 R, co e S n ~ l . 

p=0 l=-p 

Then we can write the expression on the right-hand side of (2.4) as 

co p 

(2.7) £ m ' n (dt,-p(p + n-2))v pl (t)v pl (t)dt. 

p=0 U-p ^ R 

As usually, denote by v" the Fourier transform of v, i.e. 

(2.8) v{y) = —= [ e~ i7 'v(t)dt, yeR. 

y2n Jr 

By the Plancherel's identity (2.7) is equal to 

oo P p 

V V r^Cry, -/» (p + n - 2)) |v(y)| 2 dy 
p=0 u- p 

co p 

(2.9) =TjTj I %e £m ' n ^' -P(P + n ~ 2 )) I^V(r)| 2 dy. 

p =o u-p ^ R 



Indeed, the imaginary part of the polynomial £, m ' n (iy, —p(p + n — 2)) is odd, while v p i is even for 
any real-valued function v p i, hence the imaginary part of the first integral above is equal to zero. 

From now we shall carry out all calculations for n = 3 (Steps II- VII) and then show that the 
general case can be reduced to the three-dimensional one (Step VIII). For brevity, X m ' 3 will be 
denoted £ m . 

Step II. We claim that X m , m e N, satisfy the following relations: 

(2.10) £ m (iy,-p(p+l)) = a p + ib p , 

where 



m-l 



(2.11) a = f[(y 2 +k 2 ), 

k=0 

m-2 

(2.12) ai = (y 2 +m 2 + l)Y\(y 2 + k 2 ), 



k=0 



(2p - l)a p _i + (y 2 + (p - 1 + mf) a p - 2 

(2.13) a p = A- * , p>2, 

y z + (m - 



REGULARITY OF SOLUTIONS TO THE POLYHARMONIC EQUATION 



and 

(2.14) b = ymY\(y 2 +k 2 ), 



m-l 



k=\ 

m-2 



(2.15) b x = ym(y 2 + m 2 - 1) Fl (y 2 + k 2 ) , 



k=i 

-(2p - l)Vi + (y 2 + (p - 1 + m) 2 ) b p . 2 

(2.16) ^ = jj- - , p>2. 

y z + (m - p) z 

When p = m and y = 0, the expressions in (2.13) and (2.16) should be understood in the sense of 
the corresponding limits. Formulas (2.1 1)— (2.16) show that each of the polynomials a m _i, a m , b m -\, 
b m contains the factor y 2 , hence the aforementioned limits are finite. Below we shall not consider 
separately the case p = m, y = 0, as all arguments can be justified for this case by a simple limiting 
procedure. 

Let us denote £ m (iy, -p (p + 1)) by £ p throughout this argument. According to (2.5) 

m-l 

-C P = (-l) m Y]((-iy + m-l-2M-i 7 + m-2j)-p(p+l)) 

7=0 
m-l 

(2.17) = (-irY](-iy + m-l-2j-p)(-iy + m-2j + p). 

7=0 

It remains to show that the polynomials above satisfy (2.10)— (2.16). First, 

m— 1 m 

Xo = (-D m [~[ ( -^ + m - 1 - 2/X-*'r + ™~ 2j) = (-1)'" n { ~ iy + k) 

7=0 k=-m+l 
m- 1 m- 1 



k=\ k=\ 



(2.18) = {-l) m (-iy)(m - iy) ]~[(-/y + k){-iy - k) = (y 2 + iym) ]~[(y 2 + k 2 ) 

and 

Xi = (-l) m Y\(-iy + m-2-2j)(-iy + m-2j+l) 



7=0 

m-2 



(-1)" 1 (m + 1 - z'y)(m - 1 - iy)(-iy)(-m - iy) ]~~[(-z'y + k)(-iy - k) 

k=\ 

m-2 

TO O 1 I 1 O 

(m - 2mzy - y - l)(miy — y ) (y + /c ) 



/t=i 
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m-2 



(2.19) = (y 2 (y 2 + m 2 + I) + miy(y 2 + m 2 - 1))Y\(7 2 + k 1 ), 



k=l 



as desired. As for p > 2, we have to show that 



r = 

y 2 + (m - p) 2 



+ ( l l _ p)2 (( 2 P ~ l )( a P~i ~ ib P-i) + (r 2 + (P - 1 + m) 2 ) (a P -i + ibp-i)) 

(2 - 20) = y 2 + (L- P ) 2 ^ 2p ~ 1)Z ^ 1 + (y 2 + (p~ 1+ -) 2 ) ^4 



It follows from (2.17) that 



m-l 



Lp-i = (-l) m Y\(iy + m-2j-p)(iy + m-2j + p-l) 

7=0 
m-l 

= (-l) m flO'r + m - 2(m - 1 - j) - j p)(/y + m - 2(m - 1 - j) + p - 1) 



7=0 
m-l 



(2.21) = (-l)»]~[( m -27-/ r - j p-l)(m-2j-/ r + j p-2), 

7=0 



and 



m-l 



£ p _ 2 = (-l)-]~[(-/ r + m-2j-/ J +l)(-/y + m-2j + / J -2) 



7=0 

m-l m-2 



(2.22) = (-1)'" f~[(-*'y + m- 2j' + j p-2) FlHy + m-2k- p-1). 

Hence, 



;=0 fc=-l 



y 2 + ( l_ p)2 (^ " + (r 2 + (p - i + -) 2 ) X,_ 2 ) 



m-l m-2 

y 2 + (m - p) 



111— l III — A. 

— — - I (-l) m Y\(-iy + m-2j + p-2)Y](-iy + m-2k-p-l)X 



j=Q k=0 

x((2p - l)(-m + 1 - iy - p) + (y 2 + (p - 1 + m) 2 ) (m+l-iy- p)) 

= Ti TV-*? + m-2k + p) ]l(-iy + m -2k-p-\)x 

y 2 + (m- p) 2 11 1 = 1 

(2.23) x(-m + 1 - iy- p)((2p - 1) - (m - 1 - /y + p)(m + 1 - iy - /?)). 
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The latter expression is equal to 

j, i H'r -m + p) 1 x 

p y 2 + (m - p) 2 {-iy + m + p) (—m + 1 - iy - p) 

(2.24) x(-m + 1 - iy - p)(-l)(-iy + m + p)(-iy + m- p) = £ p . 

This finishes the proof of (2.10)— (2.16), in particular, we have 

(2.25) %e£ m (iy, -p (p + 1)) = a p for every peNU {0}, 
and hence 

(-A) m u(x) u(x)\x\- 1 dx=yy a p (y) ^(yf dy. 

P^O U^p Jr 

Step III. In order to prove (2.2) for n = 3 one has to show that 

m m- 1 

(2.27) a q (y) > cj^y 2k + C Y\(g(G +l)-p(p+ 1)) for every 9 eNU {0}, y e R, 

k=l p=0 

with some constant C > independent of q and y. Indeed, if (2.27) is satisfied, then 

XX I a <?(?) |iv(y)| d y 

q=0 l=-q ^ R 

00 9 C ( m m ~ ' 

- c Z E I Z N* v ^ (r) l 2 + ^ FK^ + 1) - /> (p + DK/(r) 

? =0 /=- 9 ^ R U=l p=0 

oo 9 „ ( m m—\ s 

=c zz I z(^ v * z( °) : !+ Mon(^ +i )-p^ +i ^K'« * 

4=0 ^ R U=l p=0 , 

n' m m-1 N 

^ ^(^v(?,o;)) 2 + v(?,a>)]~[(-^-p(p + l))va,a;) dt 
2 u=i p=o , 

which combined with (2.26) gives (2.2) for n = 3. Let us now concentrate on (2.27). 

Step IV. The formulas (2. 1 1)— (2. 13) show that all a q are nonnegative. Moreover, since 

(2.29) (q - 1 + m) 2 > (m - qf for every q e N and m e N, 
we have a q > a q -2 for all q > 2, and therefore, 

m— 1 

(2.30) a q {y) > min{a (y), «i(y)} > f~[(y 2 + £ 2 ) for every q e N U {0}, y e R. 

£=0 



(2.28) 
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Hence, 

m 

(2.31) a q (y) > ^ y 2k for every e N U {0}, yet. 

k=\ 

Step V. As for the second term on the right hand side of (2.27), it is clear that 



m— 1 



(2.32) a q {y) > Y[{q(q + l)-p(p+ 1)) for q = 0, 1, m - 1, y 6 R, 

P =o 

since for such q the product in (2.32) is equal to 0. Now assume that q > m. If |y| > q + 1, then by 
(2.30) 

m— 1 m— 1 

(2.33) a,(y) > Y\((q + l) 2 + fc 2 ) > Y\(q(q +l)-p(p + 1)), 

/t=0 p=0 

and it remains to consider \y\ < q + 1, q > m. 
Step VI. By definition 

m-l 

**(y) = (-l) m ^e]~[(-/y + m-l-2j-^)(-/y + m-2j + ^) 

m-l 

(2.34) = ^Re ]~[(y 2 + (-m+ 1 +2j + q)(m -2j + q) + iy(2m- 4j - 1)). 

7=0 

Since q > m, we have (-m + l+ 2/ + (7)>0 and (m - 2j + q) > 0. Therefore, 

m-l 

fl«(7) = |~[(y 2 + (-m + 1 + 2 j + ^)(m - 2 j + 9 )) 

7=0 

m-l 



(2m -47-1) 
- 1 + 2j + q)(n, 

Let us show that there exists C(m) > m such that for every q > C{m) and |y| < q + 1 



(2.35) x *e n(l + ^ m ~* J ~ i; )■ 

11V y 2 + (-m + 1 + 2j + q)(m - 2j + q)> 



(2.36) *. fj(l + -y y2 - ( _ m ^ m + yi'Z -2j + q ^ \ 

A general idea is that the expression above can be written as 1 + %eR, where \R\ is sufficiently 
small when q is large. Indeed, R is a sum of at most 2 m - 1 terms, each of those being a product of 
(at least one, at most m) elements 

{2m- Aj- 1) 

(2.37) c: := iy— , 7 = 0, ...,m-l. 

7 y 2 + (-m+ 1 +2; + ^)(m-2j' + ^) J 
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If q > 2m, then |2 j - m\< q/2 and 
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2m 



< 



8m 



(2.38) \cj\<(q+l)- 

(g/2 + l)(<?/2) q 
Moreover, if q > 8m, the expression above is less than 1. Hence, 

8m 

(2.39) \%eR\ < \R\ < (2 m - 1) — . 

q 

Set now C(m) := 2 m+4 m. Then \%,eR\ < 1/2 for every q > C(m), and hence (2.36) holds for such 
q. Combining this with (2.35), we deduce that 

m-1 



1 I j 

z 4 (7) > - n(r 2 + (-« + 1 + 2j + ^)(m - 2j + ^r)) 

7=0 

1 m ' 1 /f 2^ \" 

> -n(-m + l + 2i + 9)(m-2i + 9 )>-|^J 



(2.40) 



> 



2 2 



| m- 1 

— n(^ +i )-pcp+D) 



p=0 



whenever |y] < q + 1, g > 2 m+4 m. 

Step VII. Finally, we claim that there exists Dim) > such that 



m-1 



(2.41) a q (y) > Dim) + 1) - p(p + 1)) if m < q < 2 m+4 m and |y| < q + 1. 
Since a g > a g _2 for g > 2, it is enough to show that for some D > 

(2.42) min{a m (y), a m +i(y)} > A lyl < q + 1, 
as (2.42) implies (2.41) with 



/ m-1 



(2.43) Dim) = D 



|~[(2" ,+4 m(2 m+4 m + 1) - p (p + 1)) 



Vp=0 



However, the formula (2.17) implies 



m— 1 



m-1 



(2.44) a m (0) = ]~[(l+2j)(2m-2j)>2 and a m+1 (0) = ]~[(2 + 2j)(2m + 1 - 2j) > 2, 
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so that for some e = So(m) > we have min{a m (y), a m+ i(y)} > 1 when |y| < s . On the other 
hand, by (2.30) 

m— 1 

(2.45) a q {y) > J^(eg + k 2 ) for every qeNU {0}, |y| > s . 

k=0 

Hence, (2.42) holds with D given by the minimum of 1 and the right hand side of (2.45). This 
concludes the proof of (2.41) and the proof of the Theorem for the case n = 3. 

Step VIII. Let us now consider any odd n > 3. Similarly to (2.27)-(2.28) in order to prove (2.2) 
one has to show that 

m m ~ §+2 

(2.46) %e£ m ' n {iy,-q{q + n-2)) > cJ^y^ + C f~| (q(q + n - 2) - p(p + n - 2)) 

fc=l 0—2 + 1 

P 2 + 2 

for every q e N U {0}, yeR. However, 
£ m >'\iy, -q(q + n-2)) 

n \ n 3 

= (-ir n( _/y +m_ i + 2 _2j_ ^( _zy + m + ^ ~ 2 ~ 2j + ^ 

jji-i 

= (-1)'" f^Hr + m - 1 - 2j - p)(-/ r + m - 2j + 

7=0 

(2.47) =£ m -\iy,-p(p + l)), 
for /? = q + | - |. Hence, 

£ m ' n (/y, (9 + » - 2)) = a q+ n_3 (y) 

£=1 p =o v ' 

-ZAcn((,^-^ t =-i)-(^-i)(^-i)) 

2 T 2 

(2.48) =C^y 2i + C p[ (4(4 + n- 2)- ^ + n - 2)) 
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where we used (2.27) and substitution s := p - | + |. As in (2.27) the constant C > above is 
independent of q and y. Now (2.48) leads to (2.46) and finishes the argument. □ 

3. Preservation of positivity for solutions of differential operators 

Proposition 3.1. 1 Assume that m e N and n e [3, 2m + l]nNis odd. Let us define h as the unique 
solution of 

(3.1) £ m > n (-d t ,0) h = 6, 

which is bounded and vanishes at +oo. Then, 

(3.2) £ m - n (-d t , -p{p +n-2))h>0, 
for all p e [0, m - n/2 + 1 /2]. 

Let us start with the operators of the second order. 

Let a, b, c, and d four real numbers such that b < < a and d < < c and let / be a (tempered) 
distribution on R. We define h as a solution of 

(3.3) (d t -a)(d t -b)h = f, 

and let the distribution g be given by 

(3.4) g:=(d t -c)(d t -d)h 

We want to investigate under what condition on a, b, c, and d the positivity of / entails the positivity 
ofg. 

First, one has to pick a suitable solution (clearly, by subtracting some solution of a homogeneous 
equation one can always destroy positivity). Also, / has to be reasonable in order for a solution to 
exist. 

To keep the discussion at full generality and, at the same time, ensure existence and uniqueness, 
in the preliminary results of this section we simply explicitly assume that the solution is obtained 
by convolution with the bounded fundamental solution vanishing at +oo (and, in particular, that / 
is such that the convolution integral converges), and later we will check that this hypothesis fits the 
argument proving Proposition 3.1. 

Also, we record a few conventions. As usually, a distribution / has the singular support at 
if / 6 C°°(R \ {0}). Furthermore, we say that such an / has an exponential decay at +oo if 
1/(01 ^ Ce~ at , t > 0, for some a,C > 0, and similarly, / has an exponential decay at -oo if 
1/(01 ^ CeP\ t < for some ft, C > 0. Within this section all distributions of interest are bounded 
in a punctured neighborhood of and thus, there is no loss of generality in considering all t > 
(or, respectively, t < 0) in the definitions above. 

Proposition 3.2. Fix a,b,c,d e R such that b < < a and d < < c. Let f and h be two 

distributions on R such that h is a convolution of f with the bounded fundamental solution of the 
operator (d t — a) (d t — b), which vanishes at +oo. In particular, 

(3.5) (d t -a)(d t -b)h = f. 



The idea of the proof has been suggested by Marcel Filoche. 
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Assume that f is positive. If [b, a] c [d, c], then 

(3.6) g:=(d t -c)(d t -d)h 

is also a positive distribution. 

Furthermore, if f has the singular support at and exponential decay at +oo and —oo, then 
g also has the singular support at and exponential decay at +oo and -oo provided that either 
a,c ± or a = c = 0. 

Proof. Let us start with the case / = 6. Problem (3.5) therefore writes 

(3.7) (d t -a)(d t -b)h = S. 

The solution of (3.7) which is bounded and vanishes at +oo is given by 

(e at if f < 1 

(3.8) h (t) = K\ bt with k=- . 

\e if f > b-a 

This solution is unique because any other bounded solution vanishing at +oo, h u would be such 
that h\ - h solves the homogeneous problem corresponding to (3.5). Thus, hi - h would have an 
exponential growth at -oo and/or violate the assumption of vanishing at +oo. 
Applying the operator in (3.6) to this solution h yields 



(3.9) g (t) = S + 



\ k (a - c)(a - d) e at if t < 0, 

k (b - c)(b - d) e ht if f > 0, 



i.e., go is the sum of the Dirac delta function, coming from the second derivative appearing in (3.6), 
and a piecewise continuous function. The constant k being always negative, it follows immediately 
that go will be positive if and only if (a - c)(a - d) and (b - c)(b - d) are both non-positive, i.e., if 
both a and b belong to the interval [d, c]. On the contrary, due to the Dirac distribution present in 
go, if either a or b is outside the interval [d, c] then g has no specific sign. 

We now turn back to the initial problem (3.5). As per our assumptions, we consider the solution 
of this problem obtained from a convolution product of / with the fundamental solution h defined 
in (3.8), that is, 

(3.10) hit) = ( f(u)h (t-u)du. 
Applying the operator in (3.6) to this h yields 

(3.11) g(t)= I f(u)g Q (t-u)du 

Jr 

We see here that the positivity of g for any positive / follows from the positivity of g , hence a and 
b belonging to the interval [d, c]. Given that b < < a, this amounts to [b, a] c [d, c]. 

Furthermore, if / has an exponential decay at both +oo and -oo and either a, c^0ora = c = 
then, due to formula (3.9), go also has an exponential decay at both +oo and -oo, and hence, so 
does their convolution in (3. 1 1), as desired. □ 
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Corollary 3.3. Let P and Q be two polynomials of order 2m, with all their roots real and their 
higher order coefficient equal to 1. Therefore P and Q write 

m- 1 m- 1 

(3.12) P(x) = Y\(x-a i )(x-b i ) and Q(x) = ]^(x - c t )(x - d t ). 

(=0 i=0 

Let us assume further that both P and Q have m non-negative roots, noted respectively (a,) and 
(c,) (0 < i < m - \), and m negative roots, noted respectively (bj) and (dj) (I < j < m), such that 
Vz e [0,m - 1], [bi,cij\ c [d { , c,]. Moreover, either none of the roots of both polynomials is zero, or 
both P and Q have exactly one zero root a = c = 0. 

Let f be a positive distribution with the singular support at and exponential decay at +oo and 
-oo, and h be the unique bounded solution of 

(3.13) P(d t )h = f, 
which vanishes at +oo. Then 

(3.14) g:=Q(d t )K 

is also a positive distribution with the singular support at and exponential decay at +oo and -oo. 

Remark. The gist of this Proposition, as well as Proposition 3.2, is preservation of the positivity 
under appropriate change of a differential operator, or more precisely, under the change that widens 
the distances between the roots. The exponential decay is not a necessary condition, but rather a 
convenient restriction that ensures the choice of a suitable solution at each step. 

Proof. First of all, let us mention that there exists a solution h of (3.13) which is bounded and 
vanishes at +oo. It is the solution obtained by the convolution of / with the fundamental solution 
of the operator P(d t ) which is bounded and vanishes at +oo. 

The Corollary will be obtained by successively applying Proposition 3.2 for all pairs [Z?„ a,-]. The 
function f Q is defined as / = /, and fi+\ is deduced from fi by 

(3.15) f M -^(dt-cdidt-ddhi, 
where h, is the solution of 

(3.16) (d t -a i )(d t -b i )h i =f i , 

obtained by convolution with the bounded fundamental solution vanishing at +oo, following Propo- 
sition 3.2. Then it is immediate to see that, by Proposition 3.2, all f are positive distributions with 
the singular support at and exponential decay at +oo and -oo. 
We claim that 

i- 1 m- 1 

(3.17) hi = Y\(d t - cj)(d t - dj) Y] (dt ~ aj){d t - bj)h, i = 1, m-1, 

y=o j=i+i 

and 

(' m- 1 

(3.18) f M = Y\(d t - Cj )(d t - dj) Y] (d t ~ aj)(d t - bj)h, i = 0, ...,m - 1, 

7=0 j=i+\ 
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where in the definitions of h , /z m _i, f m one of the products is void. 
First of all, letting f = f, we see that 

m— 1 

(3.19) h' :=Y\(d t -a } Xd t -bj)h, i = \,...,m-2, 

j=i 

satisfies equation (3.16). The only question is whether it is the same solution as ho- However, 
much as for the operators of the second order in Proposition 3.2, h' Q is bounded and vanishes at 
+oo for a particular h that we are working with (see the beginning of the proof). Clearly, the same 
is true for h obtained from (3.16) by convolution of / with the bounded fundamental solution 
of (d t - cio) (d t - bo) which vanishes at +oo. However, then the difference h - h' Q must solve the 
homogeneous equation for the operator (d, - a ) (d, - b ) and at the same time, h - h' Q must be 
bounded and vanish at +oo, hence, ho-h' = 0. Note, further, that f\ can now be defined equivalently 
by (3.18) and (3.15) and, moreover, f x has an exponential decay at both +oo and -<x>. The latter 
follows from (3. 15)— (3. 16) and the last statement of Proposition 3.2. 

Following the same lines, we verify (3. 17)— (3. 18) identifying them with functions in (3.15)— 
(3.16) as the unique solutions of the corresponding differential equations satisfying aforementioned 
decay properties. 

As we already pointed out, by virtue of Proposition 3.2, all /; are positive distributions with the 
singular support at and exponential decay at +oo and -oo. It is now clear from (3.18) that f m = g. 
This finishes the proof of the Corollary. □ 

Proof of Proposition 3.1. We now focus on the following differential operator 

£ m ' n (-d t ,-p(p + n-2)):= 

m— 1 



m - - + - - 2 j - p\\d } + m + - - ~- - 2 j + p 



m— i i 

(3.20) (-D m ]~[ [dt + 

7=0 ^ 

Let us introduce the numbers 

(3.21) Cj = 2j-(m- r ^-\, 0<j<m-\. 



Note that 

n 3 n 1 

m + - - - - 2j + p = -m + - + - + 2(m - 1 - j) + p. 



The operator X" 1 '" thus writes 



m—l 



(3.22) (-d t , -p(p +n-2)) = (-1)'" ]~[ (d t - c } - p) (d t + c j + l+p). 

7=0 

The roots of the corresponding polynomial are Cj + p and —cj - 1 - p, j = 0, m-1. 

n- 1 

Consider an integer p such that < p < m — =: k and let us examine the following set of 

roots: 

(3.23) S = {cj}U{- Cj -l} ={-k,-k-2,...,2(m-l)-k} 
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(3.24) 



U{-2(m- \) + k-\,-2(m- l) + k+ 1, ...,*-!}. 



This set contains 2m elements. The largest of the {c 7 } is larger than the largest of the {—cj - 1}: 

(3.25) 2(m-\)-k>k-\ because 2(m - 1) > Ik > 2k - 1. 
Similarly, the smallest of the {-cj - 1} is smaller than the smallest of the {cj}: 

(3.26) -2(m- l) + k- 1 < -k because 2(m - 1) > 2k > 2k - 1. 

Moreover, the {cj} and the {-cj - 1} have different parities. So when these two subsets of roots 
overlap, the overlap covers all integers, odd and even. Hence, one can rewrite the set S as the union 
of 3 sets, the center interval and two "lateral" sets, these two being respectively composed of only 
negative and positive numbers: 



(3.29) {k + 2,k + 4,...,2(m-l)-k}. 

The first and the third set might possibly be empty. 

Next, the set of roots S p = {cj + p\ U {-cj -p—1) can be described in similar terms replacing k 
by k - p: 

(3.30) S p ={-2(m - 1) + k - p - 1, -k + p - 5, -k + p - 3}U 

(3.31) ([-k + p - 1,/fc-p] nZ)u 

(3.32) {k - p + 2, k - p + 4, 2(m -l)-k + p\. 

To compare S and S p , one rewrites So in the following way 

(3.33) S ={-2(m -l) + k-l, -k -5,-k-3} U ([-it - 1, -Jfc + p - 2] n Z) 

(3.34) ([-£ + /?- l,t-p]nZ)u 

(3.35) ([k-p+ U]nZ) U {Jk + 2,Jk + 4,...,2(m- 1)-^}. 

Now the center interval is identical in both sets and composed of consecutive integers, k - p + 1 
negative numbers and k - p + 1 non negative numbers. In particular, this interval always contains 
{-1,0}. Also, compared one-by-one, the roots in the first subset of S p (3.30), are smaller than the 
roots in the first subset of 5 (3.33), and similarly, the roots in the last subset of S p (3.32), are 
larger than the roots in the last subset of So (3.35). 

We now can name the sequence of roots of S (in the order they appear in (3.33)— (3.35)) 

(3.36) b m -\ < b m -2 < ... < bo = -1 < do = < ai < ... < a m , 

and, respectively, the roots of S p (in the order of appearance in (3.30)-(3.32)) 

(3.37) d m -i < d m -2 < ■■■ < d Q = -1 < cq = < c\ < ... < c m . 



(3.27) S ={-2(m- !) + &-!, 



-k - 5, -k - 3} U 



(3.28) 
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Then \b h a,] c [d h c,] for all i = 0, m - 1, and moreover, a = c = are the only zero roots. 
Thus, we are in the setting of Corollary 3.3, which finishes the argument. (Note that when m is 
odd, we apply Corollary 3.3 with -h in place of h). □ 

4. Integral inequalities and global estimate: the case of odd dimension. Part II: weight g 

Incorporating the results of Section 3, in this section we introduce a new weight function g. 
Together with the inequalities in Section 2, this yields further improvement of the key integral esti- 
mates and ultimately, improved pointwise estimates on the solution of the polyharmonic equation. 

We start with the following auxiliary result which provides an explicit formula for the solution 
of (3.1). 

Lemma 4.1. Assume that m e N and n e [3, 2m + 1] n N is odd. Consider the equation 
(4.1) £ m '"(-d t ,0)h = 6, 

where 6 stands for the Dirac delta function. A unique solution to (4.1) which is bounded and 
vanishes at +oo has a form 

^ Z'jLi iij t < o. 

Here aj > 0, j = 1,2, m, J3j > Ofor j = 2, m and J3[ = are such that 



(4.2) hit) 



Ui n 1 \tn-l i n 1 \m—\ 

= H + 2 " 2 + 2 4o U ("2 " 2 + " " lj U ' 

and with the notation 

(4.4) y = (-<*!,. ..,-a m ,p U -,J3m), 2= (v u v m , -flu -/x m ) 



the coefficients Vj, fij e R satisfy 

(4.5) K i = (-ir +l ([\<y J -ydy 1 . 



Proof. Recall the representation formula (3.22) and let p = 0. Similarly to the analysis in the proof 
of Proposition 3.1, one can see that the characteristic polynomial of the differential equation (4.1) 
has 2m distinct roots given by 

n 1 

2 ~ 2 

and 



(4.6) cj := -m + - - - +2j, j = 0, 1, ...,m - 1, 



(4.7) ----+m-2j = -l-(-m + ---+2j) = -l-Cj, j = 0, 1, ...,m - 1. 

Following (3.27)-(3.29), (3.33)— (3.35) and (3.36), we conclude that there are m distinct negative 
roots, 1 root equal to zero and m- 1 distinct positive roots for the characteristic polynomial. Hence, 
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a bounded solution vanishing at infinity must have the form (4.2). In addition, h should be such that 
that d\g are continuous for k = 0, 1, 2m - 2 and lim f _ + df n ~ l h(t) - lim^ - df n ~ l h(t) = (-1)"'. 
This gives rise to the system of equations 

m m 

(4.8) Yj v j<~ a jf ~ Z^Oe/ = 0, k = 0,l, 2m - 2, 

;=i ;=i 

m m 

(4.9) ^ vji-ajf"- 1 - £//,(£'/)-"' : = (-l) m . 

7=1 7=1 

With the notation (4.4) and the convention 0° = 1, (4.8)-(4.9) can then be written as Ak 1 - = 
(-1)'"(0, ...,0, 1)" 1 , where A is a matrix of elements Ap = yT , i,j = 1, ...,2m. Now using the 
formula for the Vandermonde determinant, we find that the solution is k = (Kjfj™ v 

(4.10) kj = (-l) j+m Y\ ( 7j - yd (Y\(7j ~ 7d)~ l = (-l) m+1 {Y\(7j ~ yd)'' ■ 

i<k, M j, k+ j i<k )+i 

□ 

Theorem 4.2. Assume that m e N and n e [3, 2m + 1] n N is odd. Let Qbe a bounded domain in 
W, OeR n \Q,ue C™(Q.) and v = e( m -^'(u o x~ l ). Then for every £ eQandr = log l^l" 1 

(4.11) f v 2 (r,co)dco<C \ (-A) m u(x)u(x)g(log\x\'\log\^)dx, 
with 

(4.12) g(t,T) = e t (C 1 h(t-T) + C 2 ), t,reR, 

and h given by Lemma 4.1. Here C, C\, C2 are some constants depending on m and n only. 

Proof. Step I. First of all, let us pass to (t, to) coordinates. We aim to show that for some positive 
constants C and C 

\ v 2 (T,co)du) < C \ \ £ m ' n (d t ,6)v(t,u))v(t,co)h(t-T) dcodt 

(4.13) +C \ \ £ m ' n (d t ,6)v(t,co)v(t,co)dcodt. 

Jr Js--' 

Clearly, (4.13) implies (4.11). 

For future reference, observe that due to (4.1) 

(4.14) £ m ' n (-d t , 0)h(t - t) = 6(t - t), f , t e R. 

Step II. Going further, let us concentrate on the first integral on the right-hand side of (4.13): 

£ m ' n (d t , S)v(t, to) v(t, co)h(t - r) dcodt 

S n-1 
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(4.15) 



-Id 

Jr Js n - 
Jr Js"- 



£ m '"(0, 5)v(t, to) v(t, co)h(t - t) dtodt 



(£ m ' n (d t , 6) - £ m '"(0, 6)) v(t, to) v(t, to)h(t - r) dcodt. 



All terms of the operator (£ m '"(d t , 5) - £ m ' n (0, 5)) contain d k t for some k > 1, and therefore, we 
can write 

(4.16) £ m ' n (d t ,6)-£ m ' n (0,6) = d ik (-6yd*, for some d ik e R. 
Hence, 



*> 1 , i>0 
2i+k<2m 



I! 

Jr Js 



(£ m ' n (d t , 8) - £ m ' n (0, 8))v(t, to) v(t, to)h(t - r) dtodt 



(4.17) 



= V dik I I d k t Vy(t,to)Vy(t,to)h(t - t) dtodt. 
t trL Jr Js><- 1 



k>l,i>0 
2i+k<2m 



We claim that 



(4.18) f d k t vvhdt=\ f v 2 (-d,) k hdt+ V b\ f (djv) 2 djhdt, b^eR, 
Jr 2 J R J R 

for any > 1 . This can be proved by induction. For k = 1 and k = 2 we have, respectively, 

f d t vvhdt= \ f v 2 (-d t )hdt, 
Jr 2 J R 

d 2 vvhdt=^J v 2 (-d t ) 2 hdt- J(d t v) 2 hdt. 

Let us now assume that (4.18) holds for = 1,2,...,/- 1, and prove it for = I. First, 

(4.19) f d' t vvhdt = - \ d ! - 2 (d t v)(d t v)hdt- \ d l ; x vvd t hdt. 
Jr Jr Jr 

Next, using (4.18) for / - 1 and / - 2, we deduce that the expression above is equal to 



-\ f(d t v) 2 (-d t )'- 2 hdt- V b\ 2 f(ff; l v) 2 ojhdt 

1 Jr Jr 



;1,;>o 
2i+ J</-2 
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(4.20) +1- f v 2 (-d t ) l hdt + V ^J 1 f (ff t v) 2 d } t +l hdt 

2 Jr ,*-^ n Jr 



i>1.;>0 
2r+7</-l 



and (4.20) can be written in the form (4.18) for k = I. 

Then, using (4.18) for v and V^v, (4.17) leads to the representation 

J J (£ m ' n (d t , 6) - X m '"(0, 6))v(t, co) v(t, (o)h(t - t) dcodt 
= \ y d ik f f (V w v(t,oj)) 2 (-d t ) k h(t-T)dudt 

2i+k<2m 

(4.2D + y y^ff i^vfdiKt-T) d^dt 



■h+h- 



*>l,i>0 />[, >>0 
2i+k<2m 2l+j<k 



Step II. Consider I\ first. Let us decompose v into spherical harmonics. Then 
/i = r J" d ik f f (Vlv(t,co)) 2 (-d t ) k h(t-T) dcodt 



*>l,/>0 
2i+k<2m 



Y 00 P r* 

= 2 ^ dik YiYj P i (P + n-2) i v 2 pl (t)(-d t ) k h(t-r)dt 



tel.feO p=0 /=-» 
2i'+j(:<2m 



= xEZ V J' W Z d ik pXp + n-2n-d t fh(t-T)dt 

p=0Z=-p^ R fel.feO 



2/+i<2m 

00 p 

1 V 



1 /"* 

(4.22) =-VV v 2 o/ (0 (X m '"(-5,, -/7(p + n - 2)) - £" v, (0, -/>(/> + n - 2))) /*(? - r) Jr, 

2 ^0 f=p Jr 



where we employed (4.16). 
Now, 



1 00 p p 

t Z y vJXOX'^^O, -Mi? + « - 2))/i(* - r) dt 



(4.23) =^ff X m '"(0,5)va,w)va,o;)^-T)J^. 

2 Jr Js«-i 

Also, taking into account (4.14), we see that 
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1 00 p r 1 r 

(4.24) - V V vUt)£ m ' n (-d t , 0)h(t - r) dcodt = - v 2 (r, co) dco. 

2 t^n t~„ JM. * Js"-> 



p=0 l=-p 



Hence, 



1 F /-» 

h = « Yj Z v l' (t) (£ m ' n( r d <> ~P(P + n-2))- £ m <"(-d t , 0)) h(t - r) dt 

p =o i=- p J R 

-- f f £ m ' n (0,6)v(t,co)v(t,co)h(t-T) dcodt 

(4.25) +i f v 2 (t,oj)Jo;. 

2 Js»-i 

Let us discuss first in rough terms the bounds on (4.25) and then present the details. The third 
term above isolates v 2 on the sphere, as desired (cf. (4.13)), and the second term will be combined 
with the first integral on the right-hand side of (4.15), and ultimately estimated by (2.2). As for 
the first term, we split it into the sum over p < m - n/2 + 1/2 and the complementary one. The 
sum over p > m — n/2 + 3/2 will be estimated by (2.2). The sum over p < m - n/2 +1/2 must be 
positive by itself. This amounts to showing that 

(4.26) (£ m ' n (-d t ,-p(p + n-2))-£ m ' n (-d t ,0))h(t)>0, for all t±Q, 0<p<m-^ + ^, 
or, equivalently (see (4.14)), 

(4.27) £ m ' n (-d t , -p(p + n - 2))h(t) > 0, for all t±0 and 0<p<m-^ + ^. 

This, however, is exactly the result of Proposition 3.1. 
Step III. Combining all of the above, we have 



II 



£ m ' n (d t , d)v(t, co) v(t, co)h(t - t) dcodt 



= ]- f v 2 (r, co) dco + ^ f f £ m ' n (0, 8)v(t, co) v(t, co)h(t - r) dcodt 
^ 00 p „ 



+ ^ > ; > ; I v 2 p/ (0 W n ' n (-d t , -p(p + n - 2)) - £ m \-d t , 0)) - r) dt 



fel,i>0 />l,j>0 
2i+k<2m 2l+j<k 



+ Z E *S I (^Kvfdihit-r) dcodt 



> ^ f v 2 (r, w) rfw + - f f X m '"(0, <S)v(J, w) v(J, o>)/z(f - r) dc^f 



+ 
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2i+k<2m 2l+j<k 

| oo p _ 

(4.28) +2 Z Z I v^(0(X M '"(-^,- J p(p + n-2))-X'"'"(-3 f ,0))^-T)^. 

p=m-§+§ '=-/> 

Note that X m "(0, (5) is exactly the operator showing up in the last expression in (2.2). Indeed, 



/ i o \ 2 2 

(4.29) r"(0,«5) = (-If [] (m Y--2j)( m +~--2i)+5J = [] (S u - pip + n - 2)) . 

Then the absolute value of the second integral on the right-hand side of (4.28) can be majorized by 
the right-hand side of (2.2). This is easy to see, for example, after the decomposition into spherical 
harmonics: 

oo q _ 

V Y £ m > n (0,-q(q + n-2))v ql (t)v ql (t)h(t-T)dt 



q=0 l=-q 



oo q „ m-jTj 

~ZZ I n (^ + n-2)-p(p + n-2))(v ql it)) 2 \h(t-r)\dt 



?=o '=-? p=-'i+i 



(4.30) 



^ZZ J P| (^ + n-2)-p(p + n-2))(v,K0) 2 ^ 

oo ^ 

= C ZZ -C m, "(0 ! -?(? + «-2))v, / (?)v fl /(?)^ 

q=0 l=-q ^ R 

since /i is bounded by a constant and the polynomial 

(4.31) P| (q(q + n-2)- p(p + n-2)) >0 for all ?6NU(0). 

p=-|+| 

Hence, (4.1 1) follows from (4.28) once we prove that 



y y \bfj\ f f {d^vfWMt-r^dojdt 



k>l.i>0 l>l.j>0 
2i+k<2m 2l+j<k 



(4.32) 



<c r (- 

Jr" 



A) m u(x)u(x)\x\ 1 
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and 



v 2 M) \U m '"(-d„ -Pip + n - 2)) - £ m -\-d t , 0)) hit - r)| dt 



(4.33) <C \ (-A) m u(x)u(x)\x\- 1 dx 
for some C > 0. 

Step IV. Since all derivatives of h are bounded (see (4.2)), 

y y ifcgi f f (^v>) 2 i^a-r)i^ 

*>l,i>0 />l,j>0 ^ K ^ 

2i+k<2m 2l+j<k 

(4.34) < C V f f (<9*V>(?, w)) 2 JwJ/ 1 . 



*>l,i>0 

i+k<m 



Furthermore, for all k, i as above 



(4.35) (^V>(?, oj)) 2 dojdt = J^J^(p(p + n- 2))' {d k t v pl {t)f dt 



p=0 I=-p 



We break down the sum above into two parts, corresponding to the cases p < m - \ + \ 
p > m - | + |, respectively. In the first case, 

m-f+| p m-f+| p 

J £(pQ7 + /z-2))' \{%v pl {t)) 2 dt<C m , n Yj Yj (dNpKO) 2 ^ 
p=o ;=-p ^ R p=o i=- p J R 

(4.36) <C m ,„ | f id*v(t, co)) 2 dcodt < C m ,„ f (-A)'V-*M*)l-*f 1 dx, 

JrJs"- 1 Jr" 



where we used (2.2) for the last inequality (note that k > 1 since we are in the range of (4.35)) 

n , 3 
2 2' 



As for p > m - f + |, we claim that 



(4.37) (p(p + n- 2)) m <C \~\ (p(p + n-2) - s(s + n- 2)), for every p > m - § + §, 

1—-2 + 2 

2 2 

where C > depends on m and only. Indeed, one can choose C such that 

< FT L_sis + n-2)\ 
11 \ pip + n-2)y 



(4.38) i< 



fj (m-f + l)(m + |-|) ) 
( m _| + |)( m + | _I) 



for every /? > m - | + |. 
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However, since i + k < m then by Young's inequality 

(p(p + n- 2))' f (d k t v pl (t)) 2 dt= f(p(p + n- 2))V 2 "| : v(r)l 2 d J 
Jr Jr 

< f((p (p + n - 2)f + r 2m ta r )| 2 dy 

Jr 

(4.39) < (p (p + n - 2)f JV P /(0) 2 dt + J^v^)) 2 dt. 



Therefore, 



co p „ 

J] Yj{ p{p + n - 2) ) ^v P i(t)) 2 dt 

l=-n ^ R 



co p m ~i + i 



^ r 2 2 /-» 

- Z Z [I (P + » - 2) - J (J + n - 2)) (v p/ (0) 2 A 

ij 3 /—_,-, >; 3 »^R 

p=m-j+j '- P .s=-7 + § 

co P p p 

(4.40) + V V I (d'"v p i(t)) 2 dt<C I (-ArMWM(x)Ur 1 ^, 

~f 3 r - 1 Jr Jr" 



by Theorem 2.1. 

Analogously, the left-hand side of (4.33) will give rise to terms 

co p p 

(4.41) Yj Z(^ + n-2) )' I v p' WA 
for < i < m, which are bounded by 

oo p ™.-\ + \ 

C Z Z 11 {p(P + n-2)-s(s + n-2)) (v pl (t)) 2 dt 

(4.42) <C | (-Af^WMWIJcr 1 ^, 

JR" 

as above. □ 
5. Integral identity and global estimate: the case of even dimension. Part I: power-logarithmic 

WEIGHT 

Theorem 5.1. Assume that m e N and n e [2, 2m] n N even. Lef Qbe a bounded domain in R", 
O e R" \ Q, u 6 C^°(fi) anJ v = e( m ~5)'(w o x -1 ). Furthermore, let R be a positive constant such 
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that the support of u is contained in B 2R , C r := \og(AR), and let if/ be a weight function such that 
either \f/{t) = C R + tfor all t € R or if/(t) = 1 for all fell Then whenever m is even, 

Xm m-k n f 2 
(-A) m u(x) u(x) «Klog \x\~ l )dx >CVV I (dfV>) tKO dweft 
ti to Jr 

(5.i) +c J~X v ri(~ <5w ~ p ^ ,+,i ~ 2 ^ 2v ^ rf ^' 

where the product is over p = -n/2 + 2, -n/2 + 4, m - n/2 - 2, m - n/2, that is, p = -n/2 + 2j 
with j = 1,2, m/2. Ifm is odd, 

m m-k 



(-A) m u(x)u(x)^(\og\x\- 1 )dx>cyy I I (d?V>) i/j{t)dcodt 
(5.2) +C f f v n(-^-p(p + n-2)) 2 (-^ + (n/2-l) 2 )v(A(0^^, 



where the product is over p = -n/2 + 3, -n/2 + 5, ...,m-n/2-2,m-n/2, that is, p = -n/2+ 1 +2j 
wzY/z j = 1, 2, (m - l)/2. In both cases C > is some constant depending on m and n only. 

Proof. Passing to the coordinates (t, co) and decomposing v into spherical harmonics, we see that 

\ (-A) m u(x)u(x)i//(log\x\- 1 )dx= \ \ £™' n (d t , 6Jv(t, co) v(t, co) i//(t) dcodt 

Jr" Jr Js"- [ 

(5.3) = E Z £ °' n(dt > ~ q{q + n ~ 2))v ^ t] V « l(t) m dt > 

.. i\ „ ^R 



q=0 l=-q 

where 

m-l 



/ \ 

Xr(^, <U = (-l) m f] (-3, + m - n/2 - 2y)(-«9, + m + n/2 - 2 - 2j) + <5j 

;=() V / 
m-l / \ 

(5.4) = \~]l-(d t -m + 2j+l) 2 + (n/2-lf -6J. 

y=o V / 

Let us denote A := (n/2 - l) - Sco- Then the expression above is equal to 

m-l 

[~[( VA - d t + m - 2j - l)( VA + d, - m + 2j + l) 

7=0 

m- 1 m- 1 

(5.5) = ]^[(VA-d, + m-2j- l)(VA + d, + m-2j- l) = + ( Va + m - 2j - l) 2 ). 
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Hence, 

m- 1 

(5.6) £^\d t , -q(q + n-2)) = Y\ {~d 2 + Sj(q) 2 ) , 
with 

2 2 

(5.7) Sj(q) 2 = (V(«/2- l) 2 + q(q + n-2) + m-2j-\) = (q + n/l + m-lj -2) , geNUjO}. 
We claim that 

(5.8) Sj(q) 2 >Cq(q + n- 2), when q e N U {0} is such that q 2j - m - n/2 + 2, 
for some C > depending on m and n only, and 

(5.9) Sj(q) = if q = 2j -m-n/2 + 2. 

Indeed, if q > m - n/2 then (5.8) holds with any C e (0, 1) satisfying 

(5.10) (1 -C)q(q + n-2) > (m - n/2)(m + n/2 - 2) Vq > m - n/2, q e N. 

Furthermore, one can directly check that &j(q) = for the values of q described in (5.9) and 
otherwise !Bj(q) 2 > 0. Hence, if we choose C to be the minimum of !Bj(q) 2 / (q(q + n - 2)) over all 
j = 0, ...,m - 1, < q < m - n/2, such that q ± 2j - m - n/2 + 2, we obtain (5.8) for the values of 
q < m - n/2. This finishes the proof of (5.8). 
Now one can show that 

f £™> n (d t , -q(q + n- 2))v g[ (t) v ql (t) ifr(t) dt 
Jr 

= X y [(-d 2 ) m - j S k M 2 --^k i (q)\i(t)v ql (t)i/ / (t)dt+ f(-d 2 ) m v ql (t)v ql (t)mdt 

S kl (q) 2 ...S kj (q) 2 (dT j v ql (t)) 2 if,(t) dt+ f {-d]Tv ql (t) v ql (t) ifr(t) dt 

Jr 

X r m- 1 
(q(q + n-2)y(# t v ql (t)fiKt)dt+ I f[ Bjiqfvfct) ifct) dt 
... .. . ^ R ;'=0 

(5.11) >CV [ (qiq + n-2)) i (^v ql (t)fi/ f (t)dt+ f ^(O, -q(q + n- 2))v 2 ql (t) iff{t) dt. 



m n 

Zip 



/>0, k>\ 

i+k<m 



Note that the second equality uses the particular form of the weight function. Indeed, when iff = 1, 
this is just integration by parts, and when iff(t) = C R + 1, t e R, we prove by induction in k 6 N U {0} 
that 

(5.12) \ {-d 2 ) k v 2 (t) (C R + t)dt= \ {d k t v) 2 {t) (C R + t) dt, keN. 

Jr Jr 
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Finally, a straightforward calculation shows that whenever m is even, 

(5.13) £':' n (0,-q(q + n-2)) = Y\(q(q + n-2)-p(p + n-2)) 2 , 9 eNU(0(, 

p 

where the product above is over p = -n/2 + 2, -n/2 + 4, m - n/2 - 2, m - n/2, and if m is odd, 

(5.14) £7(0, -q(q + n-2)) = \\ (q(q + n-2)-p(p + n- 2)f (q(q + n - 2) + (n/2 - l) 2 ) , 

p 

for every q e NUjO}, with the product above over p = -n/2 + 3, -n/2 + 5, ...,m-n/2 - 2, m -n/2. 
This finishes the argument. □ 

6. Integral identity and global estimate: the case of even dimension. Part II: weight g 

Lemma 6.1. Assume that m e N, n e [2, 2m] n N is even and m — n/2 is even. Recall that in this 
case 

m-1 



(6.1) £':< n (-d t ,0) = Yl + (m - \ - 2j) 2 ) 



and consider the equation 

m-1 

(6.2) ^](-ai} + (m-j-2j) )h = 6, 

7=0 

where 6 stands for the Dirac delta function. A unique solution to (6.2) which vanishes at +oo and 
has at most linear growth or decay at — oo has a form 

z M2)/2 v (l) e -2it + Z (m-„/2)/2 y (2) ^ + ^n/2-1 y (3) ^m-n/l^t ^ f > Q> 

^m-n/2»2 ^(1) ^ + ^(m-n/2)/2 fg2ft + ^n/2-1 ^(3) g(m -„ /2+2i> + ^ + ^ ^ f < Q 



(6.3) = 



//ere Vj , v\ \ \ i = 1, (m - n/2)/2, v\ , z = 1, ...,n/2 - 1, and , jx (5) are some real 

numbers depending on m and n only. 

Proof. The characteristic polynomial of the differential equation (6.2) has n/2- 1 single roots given 
by 

(6.4) even numbers from - m - n/2 + 2 to - m + n/2 - 2, 
then m- n/2 + 1 double roots given by 

(6.5) even numbers from -m + n/2 to m-n/2, 
(including 0), and then another n/2 - 1 single roots given by 

(6.6) even numbers from m-n/2 + 2 to m + n/2 -2. 

This determines the structure of the solution (6.3). Furthermore, the solution to (6.2) h should be 
such that that d k t h are continuous for k = 0, 1, 2m - 2, i.e., lim^ + <9f/i(f) - lim^ - d k t h(t) = for 
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k = 0,1, 2m - 2, and lim,_o+ dj m ~ l g(t) - lim^ - dj m ~ l h(t) = (-1)'". Hence, the coefficients in 
(6.3) satisfy the linear system 

(6.7) A/? x = (-l) w (0,...,0,l)\ 

where the vector on the right-hand side is of the length 2m (it has first 2m - 1 entries equal to 

and the last entry equal to 1), the vector k is also of the length 2m with the entries (in order) 
vf\ vf\ i=l, (m - n/2)/2, vf\ i = 1, n/2 - 1, and -/uf\ i = 1, (m - n/2)/2, -/uf\ 

1 = 1, n/2 - 1, -// 4) , -// 5) . Finally, A is a 2m x 2m matrix built as follows: 

A^ = (-20 7_1 , forz' = /c, ;'=1,. ..,2m, /c = 1, (m - n/2)/2, 
A 0k = 0, for Jc = (m- n/2)/2 + 1, ...,m - n/2, 

A# = (j - l)(-2z') J "" 2 , for i = k-(m-n/2)/2, j = 2,...,2m, k = (m - n/2)/2 + 1, ...,m - n/2, 
A 7 /t = (-m + n/2 - 2i) J ~ , for i = k - m + n/2, j = 1, 2m, k = m - n/2 + 1, m-1, 
A y -£ = (2iy~ , for i = k — m + \ , j = 1, 2m, k = m, (3m - n/2)/2 - 1, 
A 0j t = 0, for = (3m - n/2)/2, 2m - n/2 - 1, 

A y * = (j - l)(2zy'" 2 , for i = k-(3m- n/2)/2 + 1, j = 2, 2m, k = (3m- n/2)/2, 2m - n/2 - 1, 
A y -£ = (m - n/2 + 2z') 7_1 , for i = k - 2m + n/2 +1, 7 = 1, 2m, /c = 2m - n/2, 2m - 2, 
Ajfc = 0, for 7 = 1 and j = 3, 2m, and A 2 u = 1, when /c = 2m - 1, 
Ajk = 0, for j = 1, 2m, and Au = 1, when /c = 2m. 

It remains to show that the system (6.7) has a unique solution, i.e., that the determinant of A is 
not 0. This, however, is a standard argument in the theory of ordinary differential equations since 
determinant of A is Wronskian at for a complete system of solutions for an ODE of order 2m. 
We omit the details. □ 

Lemma 6.2. Assume that m e N, n e [2, 2m] n N is even and m - n/2 is odd. Recall that in this 
case 



where 8 stands for the Dirac delta function. A unique solution to (6.9) which vanishes at +oo and 
has at most linear growth or decay at — oo has a form 




and consider the equation 



m-1 




(6.9) 



(6.10) 
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h(t) 



^-n/2-l)/2 v (l) e -2it + 2 (»-»/2-l)/2 y (2) ^ + y (3) , > 0> 

E («-n/ 2 -l)/ 2jU (l) ^ + 2 (m-„/ 2 -l)/ 2jU (2) ^ + ^n/2 ^(3) g ( OT -„/ 2 -l +2! -), + ^ + ^ f < Q 



//ere vj !) , vf , [if 1 , jl^ , i = 1, (m - n/2 - l)/2, v| 3) ,//| 3) , i = 1, n/2, and fi (4) , ^ are some real 
numbers depending on m and n only. 

We would like to remark a curious difference in h for three considered situations: n odd; n even 
and m - n/2 even; n even and m - n/2 odd. While formula (6.3) is rather different from (4.2), 
the underlying differential equation in (6.1)-(6.2) is a quite natural analogue of (4.1). Lemma 6.2, 
however, deals with a different equation (6.8). Just as in Section 4, we will show below how one 
can build a suitable weight from the corresponding function h. 

Proof. The characteristic polynomial of the differential equation (6.9) has n/2 single roots given 
by 

(6.11) even numbers from - m - n/2 + 1 to - m + n/2 - 1, 
then m- n/2 double roots given by 

(6.12) even numbers from -m + n/2+1 to m - n/2 - 1 , 
(including 0), and then another n/2 single roots given by 

(6.13) even numbers from m-n/2 + 1 to m + n/2-1. 

In fact, these are the same collections of numbers as those appearing in (6.4)-(6.6), in the fol- 
lowing sense. Let us denote the sequence of numbers described in (6.4)-(6.6) by S e mn and the one 
described by (6.1 1)-(6.13) by S° m>n . Then 

iS e ) d is" 1 

I m '"JmeN,ne[2,2m]nN,m-n/2even I m '"JmEN,n£[2,2m]nN,m-«/2odd ' 

via the identification S° nn = S e mn+2 for m, n as above. This identification immediately leads to the 
observation that the solution to (6.9) must be given by the formula (6.3) with n replaced by n + 2. 
□ 

Theorem 6.3. Assume that m 6 N and n 6 [2, 2m] n N is even. Let Q.be a bounded domain in W, 
O e W \ Q, u 6 Cq(Q.) andv = e^ m ~'^\u o x~ l ). Let R be a positive constant such that the support 
ofu is contained in B2R. Then there exist positive constants C, C, C" , depending on m and n only, 
such that for every £ 6 B 1R and t = log we have 

(6.14) f v 2 (r,to)dto<C \ {-k) m u{x)u{x)g(\og\x\-\\og\Z\~ l )dx 

Js"- 1 Jw 

where Cr = log(47?) and g is defined by 

(6.15) g(t, r) = h{t - r) + yU (4) (C R + r) + C + C"(C R + t), 

where h and // 4) are given by (6.3) and (6.10), depending on the parity ofm - |. 
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Proof. Step I, the preliminary choice of weight. Recall the function h defined in (6.3) and (6.10) 
for the cases when m - n/2 is even and odd, respectively. Let 

(6.16) g(t,T) :=h(t-T)+M (4) (C R + T). t,T<=R. 



In particular, when m - n/2 is even, we have 



( Z 



(6.17) g(t,r) = 



(m-»/2)/2 y (l) e ~2i(t-r) 



i=l 



(m-«/2)/2 (2) 



= 1 



Vf> (t - T)e 



+z; 



rc/2-1 (3) -(jn-»/2+2i)(f-r) 



-2i(/-r) 



+ // 4 >(C« + t), 



? > T, 



2 (*-»/2)/2 ^(1) e 2i( t -T) + E (m-n/2)/2 ^(2) (f _ 



n/2-1 (3) (m-n/2+2/)(f-r) 



+ J u (4 >(C« + 0+j" p \ 



(5) 



and in the case when m - n/2 is odd, we have 



( I 



(6.18) g(*,T) = 



(m-«/2-l)/2 (1) -2/(f-r) 



v(m-/!/2-l)/2 (2) v . 

+ Z ;=i >f e-^-^-i+soa-T) + ^(4) (Ci? + T)> 



-2/(f-r) 



£ < T, 



? > T, 



y(m-n/2-l)/2 (1) 2i(f-r) , vn m ~ 
Zj ;= i A 4 ; e "T" Zj,= i 



M2-l)/2 jU (2) (? _ T)e2i(f - r) 



n/2 (3) (m-«/2-l+20(*-r) 



+ J u (4) (C R + 0+y« p \ 



(5) 



£ < T. 



For future reference, we record a few estimates. First, by our assumptions on £ and the support 
of u, the discussion will be naturally restricted to the case t,r > log(2i?) _1 . Hence, both C R + t and 
C R + r are positive. We remark that we do not claim positivity of h or that of ~g. However, 

(6.19) \g(t, t)| < C (m, n) + |// 4) | (C R + t), t,r> log(2/?r\ 
and 

(6.20) |5 / 3,t)|<C 1 (to,/i), *,t > logC^" 1 , 1 < / < 2m, 

for some constants C (m, n), Ci(m, n) > depending on m, n only. We note that d\ g can be defined 
at ? = t for all / < 2m by continuity, and for / = 2m one assumes t + r in (6.20). 

Step II, the set-up. We commence similarly to (4.15): 

f (-A)'M*M*)?(logW-\log|£rV*= f f £™' n (d t ,S)v(t,(o)v(t,a>)g(t,T)dcodt 
Jr" Jr Js" -1 

-C"(0, tf)v(f , w) v(J, wjgft r) dwdf 

(6.21) + f f (Xr(5^)-Xr(0,5))v(f,a;)va,a;)?a,T)^. 

Jr Js"- 1 

Furthermore, as in (4.21)-(4.22), 

CCCd,, 5) - (0, 5))v(f, w) v(f, oj)g(t, r) dwdr 



= 11 

Jr JS"- 



Jr Js" 
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E E C S f f (d l t Kv) 2 dig(t,T)dcodt=:I°+I° 2 

.Jto . / c /i I 



vj,(0 (X m,n (-3 ? , + n - 2)) - £ m -"(0, + n - 2)))g(t, r) A 



(6.22) 



+ 



2i+k<2m 21+ j<k 



for some constants depending on m, n only. 
Let us now denote 

0, when m - n/2 is even, 

1, when m - n/2 is odd. 



(6.23) p Q := 
Then 



7 i = EE ^(0 + n - 2)) - (0, + n - 2)) 

L p=Q h=-p Jr V 



-£™' n (-dt, -Poipo +n-2))\ g(t, r) dt 



(6.24) 



+ 



2X 



v 2 (t, w) 



where we used (6. 1)— (6.2), (6.8)-(6.9), as well as the fact that the operator £™' n (-d t , -p (p +n-2)) 
kills constants and thus, 

(6.25) £'n-d t , - Po ( Po +n- 2))g(t, r) = £™' n (-d t , -p (p + n- 2))h(t - r) = 8(t - r). 

Step III. The term in the representation of 1° in (6.24) associated to X" 1 '"(0, -p(p + n - 2)) and the 
first term on the right-hand side of (6.21) are, modulo a multiplicative constant, the same, and can 
be estimated as follows: 

f \ £'o'"(0, 6)v(t, co)v(t, (o)g(t, r) diodt 

oo p _ 

EE vj / (0£r(o,- J p(p+«-2))ga,T)A 

p=0 /=-p Jr 

-EE I v^OXHO' + « -2))^,t)|A 

p=0 Z=-p ^ R 

°° P rs 

-EE 4 W*(0> -p{p + n - 2)) (C (m, n) + |,u (4) | (C R + o) A 
p=0 /=- P ^ R 
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= f f X" v, (0, 6)v(t, co)v(t, co) (C (m, n) + \fi {4) \ (C R + tj) dcodt 

Jr Js"- 1 

(6.26) <C \ (-A) m u(x) u(x) (C (m, n) + |// 4) | (C R + log dx, 

Jr- 

where we used the fact that the polynomial £ m ' n (0, -p(p + n - 2)) > (see (5. 13)— (5. 14)) for the 
first inequality, (6.19) for the second one, and (5.1)— (5.2) with (5.13)— (5. 14) for the last inequality 
above. 

It remains to estimate 7° (see (6.22) for definition) and 

(6.27) VV f v 2 pl {t)(L m ' n (-d t , -p(p + n - 2)) - £^(-d t , -p ( Po + n- 2))) g(t, r) dt. 

For the terms in the expression above (and for it is sufficient to either prove positivity or, for 
those which are not necessarily positive, a bound on (the sum of) their absolute values in terms of 
the right-hand side of (6.14). 

Step IV, < p < m - n/2, m - n/2 even. At this stage, the discussion splits according to whether 
m - n/2 is even or odd. In fact, the two cases are, in some sense, symmetric, but tracking both 
at the same time would make the discussion too cumbersome. We start with the situation when 
m - n/2 is even, that is, we are in the regimen of Lemma 6.1, and then we will list the necessary 
modifications for m - n/2 odd. 

To outline in a few words the strategy, we are going to show that the terms in the expression 
above corresponding toO < p < m - n/2 even are positive by virtue of the theory developed in 
Section 3, while the terms corresponding to < p < m - n/2 odd, as well as the sum of those 
corresponding to p > m - n/2, are bounded by the right-hand side of (6.14). Essentially, for large 
p we use the same considerations as in Section 4, while for < p < m - n/2 odd we observe that 
the last integrals in (5.1), (5.2) do not vanish and provide the desired bounds. Finally, can be 
analyzed much as in Section 4. Now we turn to the details. 

Step IV(a), m - n/2 even, < p < m - n/2 even. In this case, we shall employ the strategy 
suggested by the results of Section 3. One would like to show that 

(6.28) C(A0)/! = <5 

£™'\-d t , -p{p + n- 2))h(t) > 0, t*0, for all < p < m - n/2 even. 

As we shall demonstrate below, £'o' n (-d t , -p(p+n-2)), viewed as a polynomial in d t , has a double 
root at zero for all < p < m - n/2 even and hence, the result of its action (as an operator) on ^is 
the same as the result of its action on h. Thus, it is indeed enough to consider h in place of ^in the 
terms corresponding to such values of p. 

Next, due to (5.6)-(5.7), the roots of £™' n (-d t , -p(p + n - 2)) as a polynomial in d, are given by 



+(p + n/2 + m-2j -2), j = 0, m - 1. 
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Clearly, the sequence of the roots is symmetric with respect to zero and moreover, the roots can be 
ordered in non-decreasing order as single roots given by 

(6.29) even numbers from -p-m-n/2 + 2 to p-m + n/2-2, 
then double roots given by 

(6.30) even numbers from p-m + n/2 to -p + m-n/2, 
(including 0), and then single roots given by 

(6.31) even numbers from -p + m-n/2 + 2 to p + m + n/2-2. 

Now, in the spirit of the Proof of Proposition 3.1, let us denotes the non-decreasing sequence of 
the roots of £o' n (-d t , -pip + n - 2)) for p = by 

(6.32) b m -\ < b m ~ 2 < •■■ < by < b = = a < a\ < ... < a m , 

and, respectively, the roots of £'"' n (-d t , -pip + n - 2)) for a given even < p < m - n/2by 

(6.33) d m -i < d m - 2 < ... < d x < d Q = = c < c x < ... < c m . 

It is not hard to see that for all even < p < m - n/2 there holds [&,-,flj] c [d t , c,] for all 
i = 0, ...,m - 1, and moreover, a = c = = b = d are the only zero roots. We are now in 
a position to apply Corollary 3.3. The only problem lies in extra zero roots which appear in the 
current scenario. To avoid it, we do the following. 

The goal is to establish (6.28). In the framework of Corollary 3.3, let us denote by P the differ- 
ential operator of the order 2(m - 1) such that d]P(d t ) = (-l) m £™'"(-d t , 0), / := 6, and let Q be the 
differential operator of the order 2(m - 1) such that d 2 t Q(d t ) = (-l) m £" J '"(-d ; , -p(p + n-2)) with 

< p < m -n/2 even. Then the roots of P and Q are, respectively, bj, at, i = 1, m - 1, and d u Ci, 

1 = 1, m - 1, as above. They satisfy [&,-, a,-] c [di, c,], i = 1, m - 1, and they are all non-zero. 
Hence, applying Corollary 3.3 with m - 1 in place of m we have 

(6.34) P(d t ) H = 6 => Q{d t ) H is a positive distribution for all < p < m - n/2 even, 

assuming that H is the (unique) solution to P(d t ) H = 8 which is bounded and vanishes at +oo. It 
remains to show that H = (-l) m d^h, where h is the unique solution to £,™' n (-d t ,0)h = 8 which 
vanishes at +oo and has at most linear growth or decay at -oo (that is, h is given by Lemma 6.1). 
However, evidently, d^[(-l) m h] is bounded and vanishes at +oo. Since a solution with such decay 
properties is unique, (6.34) indeed proves (6.28). 

Step IV(b), m- n/2 even, < p < m - n/2 odd. We first note that 

(6.35) £™"(0, -p{p +n-2))>C, for m - n/2 even and < p < m - n/2 odd, 

where C > is a strictly positive constant depending on m and n only. This can be seen immedi- 
ately from (5. 13)— (5. 14), or, alternatively, from (5.6)-(5.7). Indeed, the polynomial X" ! '"(0, -pip + 
n - 2)) is a product of squares. Thus, it is non-negative. Moreover, if m - n/2 even, -£'"'"(0, -pip + 
n - 2)) can only be zero when p is even too. Thus, one can take C as a (strictly positive) minimum 
of £o' n i0, -pip + n - 2)) over < p < m - n/2 odd, justifying (6.35). 
Thus, for any < p < m - n/2 odd 
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(6.36) 



£2»(-d t , -p(p + n - 2)) - £Z»{.-d u 0) g(t, r) 



< d + C 2 (C R + f) 



< £™' n (0, + n - 2)) (C 3 + C 4 (C« + 0) , t, t > log(2/?)~ 1 , 

where we used (6.19)-(6.20) for the first inequality above and (6.35) for the second one. The 
constants C ; , i = 1, 2, 3, 4, depend on m and n only. Note that the operator £g' n (-d t , -pip + n - 
2)) - £o'"(-d t ,0) has order strictly smaller than 2m, thus, the corresponding derivatives of ~g are 
continuous at t = r. 

Therefore, for any < p < m - n/2 odd we have 



(6.37) 



i 

Jr 



K,(t)\ -C'r(-d t , -p{p + n - 2)) - LTi-du 0) ) r) dt 



Jr 



< C | v 2 pl it)£'^ n i0, -pip + n - 2)) (C 3 + C 4 (C« + 0) * 



< C f (-A) m M (x) M (x) (C 3 + C 4 (Cr + log IjcP 1 )) dx, 

JR." 



as desired. 



Step V, < p < m - n/2, m - n/2 odd. The argument in this case is very similar to that in Step IV, 
except that the approach to even p for m - n/2 odd resembles the approach to odd p for m - n/2 
even and vice versa. 

Indeed, when m - n/2 is odd, the case ofO < p < m - n/2 even can be handled following 
verbatim the argument in Step IV(b), observing that 

(6.38) £'">"i0, -pip +n-2))>C, for m - n/2 odd and < p < m - n/2 even, 

where C > is a strictly positive constant depending on m and n only. Then, as in Step IV(b), 
when m - n/2 is odd, for any < p < m - n/2 even we have with p = 1 



(6.39) 



V0| -C'n-d t , -pip + n - 2)) - L m '\Su -PoiPo + n-2))\ git, r) dt 



< C J v^OXrCO, -Pip + n - 2)) (C 3 + C 4 (C R + 0) dt 

<C \ i-A) m uix) uix) (C 3 + C 4 (Cr + log IjcP 1 )) dx. 

Jr." 

Thus, it remains to treat < p < m - n/2 odd. The argument essentially follows Step IV(a). 
First of all, the roots of £^"i-d,, -pip + n - 2)) for any < p < m - n/2 odd, m - n/2 odd, 
are exactly described by (6.29)-(6.31), including the parity. Hence, for all such p the operator 
£o'"(-<9 f , -pip + n - 2)) kills constants and it is enough to show that with m- n/2 odd 



(6.40) J™' n i-d t ,-\(\+n-2))h = 8 
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£™>"(-d t , -pip + n- 2))h(t) > 0, t * 0, for all < p < m - n/2 odd. 



As before, the crux of the matter is the structure of the roots of £™' n (-d t , -pip + n - 2)) and 
£,'o' n i-d t , -1(1 + n - 2)) in the context of the results of Section 3. However, if we denote the set 
of the numbers described by (6.29)-(6.31) as S e mnp when m - n/2 is even and < p < m - n/2 is 
even and as S" nn p when m - n/2 is odd and < p < m - n/2 is odd, then they can be identified 
via S° nnp = S e mn+2p _ l . Hence, the results of Step IV(a) translate into the present setting and give 
exactly (6.40), as desired. 

Step VI, p > m - n/2. Our treatment of p > m - n/2 does not depend on whether m - n/2 is even 
or odd, and we record it here in full generality. 

First of all, by (5.8) we have Sjiq) 2 > C qiq + n - 2) for all q > m - n/2 and j = 0, m - 1, 
with C > depending on m and n only. Thus, 

(6.41) £^"(0, -pip + n-2))> Cp n \p + n- 2f\ for p > m - n/2, 

with C > depending on m and n only. Secondly, recalling the definition of po in (6.23), we have 

(6.42) LC™' n (-d t , -pip + n - 2)) - £™' n (-d t , -p (p +n- 2))j g(t, r) 

<iC l+ C 2 iC R + t))p m ip + n-2T, 
by (6.19)-(6.20). Combining (6.41) with (6.42) yields then for any p > m - n/2 



(6.43) 



oo p 

E E[ 

p=m-«/2+l I=-p JR 



v%(t)\ -Cn~d t , -pip + n - 2)) - £>n-d t , -p ipo + n-2))\ git, r) dt 



< C 



< 



oo P r 

Z If 

p=m-n/2+l l=-p JR 

E Z [ 

p=m-n/2+l l=-p UR 



v 2 p i(t)p m (p + n-2) m (d + C 2 iC R + t)) dt 



v 2 pl (t)£^ n (0, -pip + n - 2)) (C 3 + C 4 iC R + t)) dt 



<C \ i-A) m uix) uix) (C 3 + C 4 iC R + log IjcI" 1 )) dx. 



This finishes the argument for 1°. 

Step VII, the bound on 1°. The estimate on 1% is essentially already incorporated in (5.1), (5.2). 
Indeed, 



J° < 

'2 - 



z z |c s :| r r (d i t Kv) 2 \dig(t,r)\dt 

fei,/>o ;>i,j>o ^ R 



2i+k<2m 21+ j<k 
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< C V (#V>(f, cu)) 2 (d + C 2 (C R + 0) 




Jr" 

employing (6.19)-(6.20) for the first inequality and (5.1), (5.2) for the second one. 
This finishes the argument. 




□ 



7. POINTWISE AND LOCAL L 2 ESTIMATES FOR SOLUTIONS TO THE POLYHARMONIC EQUATION 

This section is devoted to the proof of Theorem 1.1. In addition, we will establish sharp local 
estimates for solutions in a neighborhood of a boundary point and "at infinity", that is, when 
moving away from a given Q e d£l. 

To start, we record for future reference a well-known result that follows from the energy esti- 
mates for solutions of elliptic equations. 

Lemma 7.1. Let Q. be an arbitrary domain in W, n > 2, Q e W \ Q and R > 0. Suppose 



for every p < 2R. 

The following Proposition reflects the rate of growth of solutions near a boundary point encoded 
in Theorems 4.2, 6.3, and ultimately provides a passage to Theorem 1.1. 

Proposition 7.2. Let Q be a bounded domain in R", 2 < n < 2m + 1, Q e R" \ Q, and R > 0. 



(7.1) (-A) m « = / in Q, /6C M (fl\B«(G)), u 6 W m ' 2 (£l). 
Then 



(7.2) 




Suppose 



(7.3) (-A) m u = f in Q, / e q°(0 \ fl«(0), mgHO). 



(7.4) 




/or every p < R, 



where C is a constant depending on m and n only, and A 



in - 




Moreover, for every x e B R / 4 (Q) n Q 



(7.6) 




< i < A. 
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In particular, for every bounded domain Q. c R" the solution to the boundary value problem 
(7.3) satisfies 

(7.7) \V n ~ n/2+l/2 u\ e L°°(Q) when n is odd and \V m ~ n/2 u\ e L°°(Q) when n is even. 
This Proposition will conclude the proof of Theorem 1.1. 

Proof of Proposition 7.2. Without loss of generality we can assume that Q = O. Let us approximate 
Q by a sequence of domains with smooth boundaries {£l n }™ =l satisfying 

oo 

(7.8) ^Jq„ = Q and fl„ c Q„ +1 for every nel 

n=\ 

Choose no e N such that supp /cfi„ for every n > n and denote by u n a unique solution of the 
Dirichlet problem 

(7.9) (-A) m u n =f in Q„, u n € W m ' 2 (Q„), n > n Q . 

The sequence {Wn}~ = „ converges to u in W'"' 2 (Q) (see, e.g., [39], §6.6). 
Next, take some tj e Cq(B 2 r) such that 

(7.10) < i] < 1 in B 2R , 7] = 1 in B R and |V fc ?7| < CR~ k , k < 2m. 

Also, fix t = logp" 1 and let g be the function defined in (4.12) when n is odd and by (6.15) when 
n is even. 

We observe that when n is odd, the formula (4.12) yields 

(7.11) \V k x g(\og\x\- l ,\ogp- l )\<C\x\- k - 1 , 0<k<2m, xel"\ {0}, p 6 (0, oo), 
while for n even by (6.15) and (6.19)-(6.20) 

(7.12) |g(log W- 1 , logp" 1 )] <d +C 2 (C R + log W- 1 ), 0<\x\,p<2R, 
and 

(7.13) IV^dogW^Jogp- 1 )! <C\x\~ k , \<k<2m, Q<\x\,p<2R. 

Here, as usually, we assume \x\ + p when k = 2m, and lower derivatives of g as well as g itself are 
defined at x such that \x\ + p by continuity. Hence, in particular, when n is odd, 

(7.14) IV^CogW-Mogp- 1 )^^-*" 1 , Q<k<2m, x e C RM , p < R, 
and when n is even, 

(7.15) IV^logUrMogp- 1 )! <CR~ k , Q<k<2m, xeC RM , p < R, 
since for x e Cr^r we have 

AR 

C R + log W" 1 = log(4tf) + log W" 1 = log — * C. 

\x\ 

Consider now 

(7.16) f ([(-A) m ,T]]u n (x))Ux)u n (x)g(\og W-Uogp" 1 )) dx, 
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where the brackets denote the commutator 

[(-A) m , rj]u n (x) = (-A) m {ri{x)u n {x)) - rj(xX-A) m u n (x), 

the integral in (7.16) is understood in the sense of pairing between W m ' 2 (Q. n ) and its dual. Evidently, 
the support of the integrand is a subset of supp V77 c C r: ir, and therefore, the expression in (7.16) 
is equal to 

(7.17) £ ([(-Ar,rj]u n (x^{r,(x)u n (x)g(log \x\~\\ogp' l )yx 

* c z w^-i f |v ' w " w|2 dx - -^-n r i"««i 2 ^ 

~^ K JC R2 R K JCrar 

using Cauchy-Schwarz inequality, (7.14)— (7.15), and (7.10) for the first inequality and Lemma 7.1 
for the second one. 

On the other hand, since u n is biharmonic in Q. n n B 4R and rj is supported in B 2 r, one can see that 
rj (-A) m u n = and hence the integral in (7.16) is equal to 

(7.18) J (-A) m (rj(x)u n (x)) r ] (x)u n (x)g(\og \x\~\ logp" 1 ) dx. 

To estimate it we employ Theorems 4.2 and 6.3 with u = rj u n . The results (4. 1 1) and (6.14) hold for 
such a choice of u. This can be seen directly by inspection of the argument or one can approximate 
each u„ by a sequence of C™(£1„) functions in W' n ' 2 (Q. n ) and then take a limit using that O £ Q. n . 
Then (7.18) is bounded from below by 

( 7 - 19 ) ^77 f \v(x)un(x)\ 2 do- x . 
P Js p 

Hence, for every p < R 



( 7 - 2 °) -p^l l""WI 2 ^<^£ \u n (x)\ 2 dx. 



'Crar 

Now the proof of (7.4) can be finished by taking the limit as n — » 00. 

Going further, by virtue of the local estimates for solutions of elliptic equations, the bound in 
(7.4) transforms into the uniform pointwise estimates for V A u (observe that the exponent A, as 
defined in (7.5), is an integer number both when n is odd and when n is even). 

Indeed, by interior estimates for solutions of the elliptic equations (see [4]) 

C 

d(x) n JB d{x)/2 (x) 

where d(x) denotes the distance from x to dQ.. Let xq be a point on the boundary of Q such that 
d(x) = \x - x \. Since x e B R / 4 (Q) n Q and QeR"\ Q, we have x e B R / 4 (x ), and therefore 

\-2i+2A 

\u(y)\ 2 dy, 



(7.21) \Vu(x)\ 2 < -^ji |V ! w(y)| 2 dy, < i < m, 



1 r r r d(xY 2l+2A C 

(7.22) \Vu(y)\ 2 dy< \u(y)\ 2 dy < C V ' - 

«W Jb Ax)/2 (x) a ^ x ) JB 2dU) (x ) K Jc 3 r,4Mxo) 
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as soon as i < A. Indeed, the first inequality in (7.22) follows from Lemma 7.1. Turning to the 
second one, we observe that d(x) < R/4 and therefore, 2d(x) < 3R/4. On the other hand, u is 
polyharmonic in B 4R (Q) n Q, and \Q — xq\ < \Q - x\ + \x - xo\ < R/2. Hence, u is polyharmonic 
in 5 3S (xo) fl O. and Proposition 7.2 holds with xq in place of Q, 3R/4 in place of R and p = 2d(x). 
Furthermore, Cs R /4^ R (xo) c C R / 4 ^ R (Q), and that finishes the argument for the second inequality in 

(7.22) . Clearly, d(x) < \x - Q\, so that (7.21)-(7.22) entail (7.6). □ 

At this point, we are ready to address the behavior of solutions "at infinity". 

Proposition 7.3. Let Qbe a bounded domain in W, 2 < n < 2m + 1, Q e R" \ O, r > and assume 
that 

(7.23) (-A) m u = f in Q, fe C™(B r/4 (Q) nQ), ue W m ' 2 (Q.). 
Then 

(1. 24) p 2 ^ 1 - 4 '" r \u(x)\ 2 do- x < Cr 2A+n ~ 4m f \u(x)\ 2 dx, 

Js p (Q)nn Jc r/4 . r (e)nn 

for any p > r and A given by (7.5). 
Furthermore, for any x e Q. \ B 4r (Q) 

j.2A+n-4m p 

(7.25) |VW)| 2 < C - „ l2A+2n _ 4m+2i I \u(y)\ 2 dy, < i < A. 



\x-Q\ 



Cr/4,4 r (2)nn 



Proof. Without loss of generality, one can consider Q = O. Retain the approximation of Q with 
the sequence of smooth domains £l„ satisfying (7.8) and define u„ according to (7.9). We denote 
by I the inversion x i-» y = x/\x\ 2 and by U„ the Kelvin transform of u n , 

(7.26) U n (y):=\y\ 2m - n u n (y/\y\ 2 ), y e I(& n ). 
Then 

(7.27) (-A) m U n (y) = \y\- n - 2m ((-A) m u n )(y/\y\ 2 ), 

and therefore, U„ is polyharmonic in I(Q. n ) fl B 4 / r . Moreover, U n e W m ' 2 (I(Q. n )) if and only if 
u n e W m,2 (Q. n ). Note also that Q. n is a bounded domain with O i Q. n , hence, so is I(Q. n ) and 
O i I(Q n ). 

At this point we can invoke Proposition 7.2 to show that 

(7.28) p 2A+ "- 1 f \U n (y)\ 2 dcr y <Cr 2A+ " f \U n (y)\ 2 dy, 

J Slip JCl/ r J/ r 

so that 

(7.29) p 2 ^' 1+1 - 4m f \ Un (x)\ 2 d<r x < C r 2A+n ~ 4m f \u n (x)\ 2 dx. 

JS p JCrjlj 

using the substitution (7.26) and the change of coordinates. 
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This furnishes the desired L 2 estimates. It remains to establish pointwise control. To this end, 
let us fix some x e Q. \ B Ar {Q) and recall that 

U n (x) = \x\ 2m ~ n U n (xl\x\\ 

hence, 

(7.30) |V*« B (x)| < C Yj W 2 ™'~* (V k U n )(x/\x\ 2 ). 

Therefore, combining (7.30) and Proposition 7.2 applied to the function U n , we deduce that 

(7.31) \Vu n {x)\ 2 < C m+2i \U n (z)\ 2 dz = C m+2i \u n {z)\ 2 dz, 

for all < i < A. Now one can use the limiting procedure to complete the argument. Indeed, 
since u n converges to u in W" l ' 2 (Q.) and i < A < m, the integrals on the right-hand side of (7.29) 
and (7.31) converge to the corresponding integrals with u n replaced by u. Turning to \Vu n (x)\, we 
observe that both u n and u are m-harmonic in a neighborhood of x, in particular, for sufficiently 
small d 

(7.32) \V(u n (x) - u(x))\ 2 < [ \u n (z) - u(z)\ 2 dz. 

" JB d/2 (x) 

Asn^> oo, the integral on the right-hand side of (7.32) vanishes and therefore, \Vu n (x)\ |V'm(jc)|. 
Similar considerations apply to u n (x) and the integral on the left-hand side of (7.29). □ 

8. Green's function estimates 

The present section addresses sharp pointwise estimates on Green's function of the polyhar- 
monic equation and the regular part of Green's function, that is, the difference between Green's 
function and the fundamental solution. We have discussed some of these bounds, as well as their 
applications for the solutions of the Dirichlet problem, in [23] . The estimates listed below are more 
refined compared to the statements in [23] and we present them here with full proofs. 

To start, let us recall the definition of the fundamental solution for the polyharmonic equation 
(see. e.g., [7]). A fundamental solution for the m-Laplacian is a linear combination of the char- 
acteristic singular solution (defined below) and any m-harmonic function in R". The characteristic 
singular solution is 

(8.1) C m; „|;t| 2m "'\ if n is odd, or if n is even with n > 2m + 2, 

(8.2) C m ,«l^l 2m " log \x\, if n is even with n < 2m. 

The exact expressions for constants C m „ can be found in [7], p. 8. Hereafter we will use the funda- 
mental solution given by 

( \x\ 2m - n , if wis odd, 

(8.3) Y(x) = C m>n J \x\ 2m - n log ^sp, if n is even and n < 2m, 

\ \x\ 2m ~ n , if n is even and n > 2m + 2. 
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As in the Introduction, by G we denote the Green function of the operator (-A) m (see (1.6)) and 
S stands for its regular part, that is S (x, y) = G(x, y) - Y{x - y), x,y e Q,. 

Theorem 8.1. Let Q. cl° be an arbitrary bounded domain, m e N, n e [2, 2m + 1] n N, and let A 
retain the significance of (7.5). Fix any number N > 25. Then there exists a constant C depending 
on m, n, N only such that for every x, y 6 Q the following estimates hold. 
Ifn 6 [3, 2m + 1] n N is odd then 

d(y) A ~ j 



(8.4) \V x Vp(x,y)\ < C _-;; + »-w when \^-y\^Nd(y), 0<i,j<A, 



and 

(8.5) lV;vjG(x,y)| < C ^^^ 2m+j , when \x-y\>Nd(x), 0<i,j<A. 
Next, 

(8.6) |V' v V^,v)| < C 



\ x — y\n-2m+i+j ' 

when \x-y\ < N~ l ma.x{d(x), d(y)}, and i + j > 2m - n, < i, j < m - n/2 + 1 /2, 



and 



\2m—n—i—j 



(8.7) \V x V J y G(x, y)\ < C mm{d(x), d(y)Y 

when \x - y\ < N" 1 max{d(x), d(y)}, and i + j < 2m - n, < i, j < m - n/2 + 1 /2. 

Finally, 

' j C C 

(8.8) x y G(x,y)\ < — — - _ y ^ n _ 2m+i+j * max{ j W; d(j ^ \ x _ y |}»-2,„ + « + >' 

when N~ l d(x) < \x - y\ < Nd(x) and N' 1 d(y) <\x-y\< Nd(y), < i, j < A. 

Furthermore, ifn £ [3, 2m+ 1] nN is odd, the estimates on the regular part of the Green function 
S are as follows: 

(8.9) \V[V(,S(x-y)\<j—^—- when \x - y\ > Nmm{d(x), d(y)}, < ij < A. 
Next, 

(8.10) |ViV^,,)| < ^ mx) C d(y)r ^, , 

when \x - y\ < N~* max{d(x), d(y)}, and i + j > 2m - n, < i, j < m — n/2 + 1 /2, 

and 



(8.11) \V x V J yS(x,y)\ < Cmm{d(x),d(y)} 



2m—n—i—j 
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when \x-y\ < N~ l max{d(x), d(y)}, and i + j < 2m - n, < i, j < m - n/2 + 1 /2. 

Finally, 

. . c C 

(8.12) |V x VyS(x,y)| < m[n{d ^ d(y ^\ x _ y^n-im+i+j ~ max {J(x), d(y), \x - y\Y~ 2m + i +r 

when AT 1 d(x) < \x-y\ < Nd(x) and N~ l d(y) < \x-y\ < Nd(y), < i, j < m-n/2+1/2. 
Ifn £ [2, 2m] n N is even, then (8.4)-(8.5) and (8.8) are valid with A = m — |, and 

(8.13) V^iG(x.y) < Cmm{d(x),d(y)} 2m - n -^ lc' + log min ^>^)} 

\ " \x-y\ 



(8.15) \V x Vi,S(x,y)\ < Cmin{d(x), d(y)} 2m ~ n ~ H [C + log 



when \x - y\ < N 1 max{<i(jc), d(y)} and < i, j < m - n/2. 

Furthermore, ifne [2, 2m] n N is even, the estimates on the regular part of the Green function 
S are as follows: 

(8.14) \V x ^ y S(x - y)\ < C\x - y^-H l C ' + log 

when \x-y\> Nmin{d(x), d(y)}, < i,j < m - n/2. 

Next, 

diam Q. \ 
max{d(x), d(y)} ) ' 

when \x-y\ < AT 1 max{d(x), d(y)}, 0<i,j<m-n/2. 

Finally, 

(8.16) \V x V y S(x,y)\ < Cmm{d(x),d(y),\x-y\} 2m - n - H {c + \og — dmmQ ,.^ 2m+ ,. + , | 

y \ max{d(x),d(y),\x - y\} n 2m+,+J ) 

when N~ l d(x) <\x-y\< Nd(x) and N~ l d(y) <\x-y\< Nd(y), < i, j < m - n/2. 

Before passing to the proof of the Theorem, we would like to point out that the bounds on the 
highest order derivatives highlighted in Theorem 1 .2 is a particular case of Theorem 8.1. 

Proof of Theorem 8.1. Recall the definition of the fundamental solution Y in (8.3). For any a, a 
multi-index of length less than or equal to A, the function d" Y(x - y) can be written as 

(8.17) d a y T{x -y) = P a (x - y) log + Q a (x - y). 

y \x-y\ 

Here, when the dimension is odd, P a = and Q a is a homogeneous function of order 2m — n — \a\. 
If the dimension is even (and less than or equal to 2m by the assumptions of the theorem) then P a 
and Q a are homogeneous polynomials of order 2m — n — \a\ as long as \a\ < 2m - n. In both cases, 
P a and Q a do not depend in any way on the domain Q. 
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Consider a function 77 such that 
(8.18) 7/eC"(fii /2 ) and 77=1 in B m , 
and define 

(x-y 



(8.19) ^(x,y):=a?G(x,y)-i7 



rf(y) 



P ff (x - y) log + Q a (x -y)}, v. v e £L 



For every fixed yefi the function x i-» ^(x, y) is a solution of the boundary value problem 
(8.20) (-AJ n K(x, y) = f a (x, y) in Q, /„(•, y) 6 ft (-,y) 6 W'"- 2 (ft), 

where 

|x-yl 



(8.21) / ff (x,y):=(-A,)"X(*,y) = 



(-A,r,/7 



P a {x-y)\og 



d(y) 



+ Q a (x-y)\. 



Indeed, it is not hard to see that for every a 

(8.22) f a (;y)eCZ(C d(y)/44(y)/2 (y)) and \f a (x,y)\ < CdiyT"^, x,yeQ, 
so that, in particular, /<*(•, y) 6 C™(Q) as stated in (8.20). Furthermore, by (8.21), 

(8.23) pPW;y)\\ LHn) = \\K(;y)\\ w ^) < Q\f (;y)\\ w - m ,2 {£1) , 
where W~ m,2 (Q) stands for the Banach space dual of W' n ' 2 (Q.), i.e. 



(8.24) \\f a (;y)\\ w - m ,2 {a) = 



sup 

velV'".2(n):||v|| lV ,„, 2(£J) =l 



f fa 

JQ. 



(x, y)v(x) dx. 



Recall that by Hardy's inequality 



(8.25) 

' -Q\* 
Then by (8.24) 



L 2 (n) 



< C||V m v|U n) for every v 6 W m ' 2 (Q), <2 6 



(8.26) \\f a (;y)\\ w - m ,2 in) 



sup 



< C sup 

veW<». 2 (n):||v|| li ,„,, 2(£i) =l 



• -yo\ r ' 



\ f a (x,y)v(x)dx 

\\f a (;y)\--yo\ m \\m a) 



< C sup d(y) m II V'"v|| i2(n) ||/ B (., y)lb ( c rf( , ww/2W) 

veW». 2 (0): 11^1^,2(^=1 



< CJ(y) m ||/ tt (-, y)|| L 2 (Q(> , )/4 ^ )/2M) , 



where y 6 <9Q is such that \y - y | = d(y). Therefore, by (8.22) 
(8.27) ||V^ a (-,y)|| L2(n) < Cd{y) m - M - nl2 . 
Now we split the discussion into a few cases. 
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Case I: either \x - y\ > N d(y) or \x - y\ > N d(x) for some N > 25. 

Let us first assume that \x - y\ > Nd(y), N > 25. In this case we have ^(^y) = and hence, 

(8.28) Y x K a (x, y) = V x dyG(x, y), for x, y e Q. such that \x - y\ > N d(y), < i, \a\ < A. 

As before, we denote by y some point on the boundary such that \y -y \ = d(y). Then by (8.20), 
(8.22) the function x \-> ^{x^) is m-harmonic in Q. \ B 3d(y y 2 (yo)- Hence, by Proposition 7.3 with 
r = 6d(y) 



d(y) 1A -' 



(8.29) iVX(^j)l 2 <C , UyyJ l? , +2n 4m+2 , f \K(z,y)\ 2 dz, for < / < A, 



provided \x - y \ > Ar = 24 <i(y). The latter condition is automatically satisfied for x such that 
\x-y\ > N d(y) with N > 25, which is the current assumption on x, y. The right-hand side of (8.29) 
is bounded by 

d(y) 2A + n-2 m n \K(Z,y)\ 2 

I n ) i O,, A ™ i 1." I i lO™ 

Jc 



v a! \2A+2n-4m+2i | |„ -., |2m 

d(y) 2A+n - 2m r . ^ J(y) 2i - 2 l ffi l 



(8 - 30) - c k-yoi^-^ J n i v r^fc^)i z * < c u v|2 . i+2 „_ 4m+2f . 

by Hardy's inequality and (8.27). Therefore, 

d(v) 2A ~ 2M 

(8.31) |VX(x,y)| 2 <C _ , 2/l+2n _ 4ffl+2 , , when \x-y\>Nd(y), 0<i,\a\<*. 

\x y\ 

By (8.28), the estimate (8.31) with j := \a\ implies (8.4). Also, by the symmetry of Green's 
function we automatically deduce (8.5). This furnishes the desired estimates on Green's function 
when either \x-y\ > N d(y) or |x - y\ > N d(x), both in the case when n is odd and when n is even. 

Turning to the estimates on the regular part of Green's function, we observe that, in particular, 
(8.4) and (8.5) combined give the estimate 

(8.32) \Y x Vi,G(x,y)\ 2 < |jc _ |2 f 4m+2 ,. +2; . , when \x - y\ > N mm{d(x),d(y)}, < i, j < A. 
Furthermore, if n is odd, then 

(8.33) \^T(x- y )\< lx _^_ 2m+i+ . forall x,y e Q, i,j>0, 
while if n is even, then 

diam (Q) 

\x-y\ 

for all x,y 6 Q. 

Combining this with (8.32) we deduce (8.9) and (8.14) in the cases when the dimension is odd 
and even, respectively. 

Case II: either \x-y\< N~ l d(y) or \x-y\< N~ l d(x) for some N > 25. 



(8.34) |V;V^r(jc -y)\ < d \x -y\~" +2m -^j log t \ / + C 2 \x-y\- n+2m - H if < i + j < 2m-n, 
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Assume that \x-y\ < N _1 d(y). For such x we have ^(^) = 1 and therefore 

(8.35) K a {x,y) = &* G(x,y) + P a (x - y) log - Q\x - y). 

y d(y) 

Hence, if n is odd, 

(8.36) <R a (x,y) = &* y (G(x,y)-T(x-y)), when \x - y\ < N~ l d(y), 
and if n is even, 

(8.37) K a (x,y) = d^(G(x,y)-r(x-y)) + P a (x-y)\og^^, when \x-y\ < N~ l d(y). 
By the interior estimates for solutions of elliptic equations for any i < m 

(8.38) \V x <R a {x,y)\ 2 < — ^ f \<R{z,y)\ 2 dz, 

since by (8.22) for any fixed y the function 7? ff (-,;y) is m-harmonic in B^/^iy) 3 B d( yy%(x). Next, 
since for every z e B diy) ^{x) c i^^Cy) and any y e such that \y - y \ = d(y) we have 
|z - y \ « d(y), one can bound the expression above by 

(o on\ C f \K(z,y)\ 2 C ■ ■ ||2 C 

1 ^ J(v)" +2l '- 2m J Brf(y)/4(y) |Z - V | 2m " d(jO" +2 '- 2 " ! 11 * 1 ' ■ y; " i2 (") - rf(y)2«-4m + 2* + 2M ' 

with < \a\ < A. Then, overall, 

(8.40) KKi^yf < d(y)2n C 4m+2i+2lar \x-y\<N- l d(y), 0<\a\<A. 

Having this at hand, we turn to the estimates on the regular part of Green's function, starting 
with the case when n is odd. It follows from (8.36) and (8.40) that 

(8.41) \V x V y S(x,y)\ < — 0<i<m, < j < A, \x -y\ < N' l d(y). 
and hence, by symmetry, 

(8.42) |V^5(j:,y)|< C 2w+/+ . , < i < A, < j < m, \x-y\ < N~ l d(x). 
However, we have 

(8.43) \x - y\ < N~ l d(y) => (N - l)d(y) < Nd(x) < (N + l)d(y), 

i.e. d{y) « d{x) ~ min{d(y), d(x)} ~ max{d(y), d(x)} whenever \x-y\ is less than or equal to N^diy) 
or N~ l d(x). Therefore, when the dimension is odd, we obtain (8.10) and (8.11) depending on the 
range of i, j. 

Moreover, by (8.33), for the Green function itself, we then arrive at (8.6), (8.7), once again 
depending on the range of i, j. 

Similar considerations apply to the case when the dimension is even, leading by (8.37) to the 
following results: 
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diam£^ 



(8.44) \V x V J y S (x, y)\<C diyf" 1 " 1 '^ \C' + log ■ 



d(y) 



forO < i < m, < j < m - -, \x - y\ < N l d(y), 



and 

diamfl 



(8.45) |V ! t VAS" (x, y)\ < C d(xf m - n - H I C + log 



d{x) 



forO < i < m - ^, < j <m,\x-y\ < N l d{x). 

Overall, in view of (8.43), and the fact that 2m - n - i - j > whenever < i, j < m - |, 
(8.44)-(8.45) yield (8.15). 

As for the Green's function estimates, when n is even, (8.35) and (8.40) lead to the bounds 

(8.46) \V x V y G(x,y)\ < C diyf" 1 ^ Ic + log , < i,j < m - \x -y\ < N~ l d(y), 

(8.47) \VMG(x, y)\ < C J(;c) 2m ^'" ; ' \C + log ) , < i, j < m - -, \x - y\ < N~ l d(x), 

■ \ \x-y\J 2 ' 

and therefore, one obtains (8.13). 
Finally, it remains to consider 

Case III: \x-y\ « d(y) « d(x), or more precisely, the situation when 

(8.48) N~ l d(x) < \x - y\ < Nd(x) and AT 1 J(y) < \x - y\ < Nd(y) for some N > 25. 

Fix any x,y as in (8.48). Then the mapping z G(z,y) is m-harmonic in Bd( x )/(2N)( X )- Assume 
first that n is odd. Then by the interior estimates, with x e dQ. such that \x - xq| = d(x), we have 



W^G(x,y)\ 2 < f |V?G(z,v)| 2 Jz 



"W JB d(x y (2 m(x) "W JB 2d(x) (x ) \Z~Xq\ 



by definition (8.19). Then the expression above is bounded by 

C f ,_l/vl„, « ,9 . C 



r iv? 

J Bd(x)/(2N)( X ) 



l .;lT(z-yfdz + ~ | |V^( Z ,y)N Z 

«W JB d(x v<2m(x) «W 



C C 

(8.50) < A9,i-4m-u?;-i-9UI 



J( JC )2«-4m+2,+2| a | d{x) n - 2m+li d(y) n - 



2m+2\a\ ' 



provided that < ?' < m, < \a\ < m - \ + 1 



2 2" 

The same estimate on derivatives of Green's function holds when n is even, since 
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\V x ^G(x,y)\ 2 < f | V NG(z,v)| 2 Jz<-|^- f \P a {z-y)Uz 



(8 - 5i) r ir(z - j)|2 * + «^ r 



\K(z,y)\ 2 
\z - x \ 2m 



since the absolute value of log is bounded by a constant for z,x,v as in (8.51), (8.48). This 
yields (8.8), both in the case when n is even and odd. 

Furthermore, if n is odd, the same argument (or simply an estimate on the difference between 
the Green function and S , that is, the fundamental solution) can be used to deduce 

(8.52) \VMS(x,y)\< C ( 



min{ d(x), d(y), \x - y\}n-2m+i+j max{d(x), d(y), \x - y \}"-^+i+j 

when x, y satisfy (8.48), < i, j < m - n/2 + 1 /2, 
which gives (8.12). However, if n is even, using (8.48), we are led to the bound 

iviv;s (v >i •: '' [r + i„ g " : " nli 1 



xy mm{d(x),d(y),\x-y\} n ' 2m+i+ j \ max{d(x),d(y),\x - y\}"- 2m+i+ j 

C / , diamQ 

(8-53) « — — - - .— _ r+log 



max{d{x), d(y), \x - y\}"-2m+i+j \ ° max {d( x ), d(y), \x - v |}"- 2m+, ' + ^ 
for < i, j < m - |. This establishes (8.16) and finishes the argument. □ 



9. Estimates for solutions of the Dirichlet problem 

The results of Section 7 provide certain local monotonicity estimates for solutions of the Dirich- 
let problem. They can, in principle, be translated into the bounds on solution u in terms of data / 
when / is a suitably supported C™ function. However, the pointwise control of the Green function 
that we obtained in Section 8 allows us to estimate the solutions of the Dirichlet problem for a 
wide, in some sense, optimal, class of data. 

Proposition 9.1. Let Q, c W be an arbitrary bounded domain, m 6 N, n 6 [2, 2m + 1] n N, and let 

A retain the significance of (7.5). Consider the boundary value problem 

(9.1) (-A) m u = ^ c a d a f a , u 6 W m \Q.). 

\a\<A 

Then the solution satisfies the following estimates. 
Ifn 6 [3, 2m + 1] Pi N is odd, then 



(9.2) 




\f a (y)\ dy, xeQ, 
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whenever the integrals on the right-hand side of (9.2) are finite. In particular, 



(9.3) W^MV-m < C m ,„, n £ l|d(O m " H " W / a |lz, ( n), P > 



W<m-|+5 



(9.4) |V m -5 M (x)| < C m ,„ >; I J(y)'"~3-I ff| log) 1 + jjj^l 1/aOOI dy, x e Q, 



provided that the norms on the right-hand side of (9.3) are finite. 
Ifn 6 [2, 2m] Pi N w even, then 

%n 2_j \ ' 

M<m-§ ^ n 

whenever the integrals on the right-hand side of (9.4) are finite. In particular, 

(9.5) l|V"^ M || L „ (n) < C m , n , n J] IW)^ Hol /«ll^(n). P > L 

M<m-§ 

provided that the norms on the right-hand side of (9.5) are finite. 

The constants C m> „ above depend on m and n only, while the constants denoted by C WiBj q depend 
on m, n, and the diameter of the domain Q. 

Proof. First of all, the integral representation formula 

(9.6) «(*) = f G(x,y)^c a <9 a / ff (y)</y, xeQ, 

follows directly from the definition of Green's function. It implies that 

(9.7) V A u(x) = ]Tc a (-\) M f V*d°G(x,y)f a (y)dy x e ft. 

|or|</l ^ n 

Let us now assume that n e [3, 2m + 1] n N is odd. We claim that in this case 

(9.8) V™ 2+2 d a y G(x,y) <C-^f — , for all x, y e ft, 0<]a]<m-- + -. 

y \x-y\ 2 2 

Indeed, when \a\ = m-| + ^, (9.8) is equivalent to (1.8). Thus, it remains consider \a\ < m—% — \. 
Now we split the cases, essentially according to Theorem 8.1. When \x - y\ > Nd(y), for some 
N > 25, (8.4) directly yields (9.8). If x, y e ft are such that \x - y\ < N~ Y d(y) then (9.8) follows 
from (8.7). Finally, if N~ l d(y) < \x - y\ < Nd(y), we observe that d(x) < \x - y\ + d(y) < 
(1 + N)\x - y\. The latter estimate shows that (8.5)— (8.8) imply the bound on the left-hand side of 

(9.8) by C \x - y\ m -i- x 2-\ a \ and hence, (9.8), since here d(y) * \x - y\. 

Now (9.8) allows us to deduce (9.2), and (9.3) follows from it via the mapping properties of the 
Riesz potential. 

Let us now turn to the case when n e [2, 2m] n N is even. Then we have to show that 

(9.9) VT l dPG{x,y) < C d(y) m ~^ lal logfl + -^-|, for all x,y eQ, < \a\ < m- -. 

y ' \ \x-y\J " 2 

Once again, we split to cases according to Theorem 8.1. When \x—y\ > Nd(y), for some N > 25, 
then 1 + ^QL is bounded from below and above by a positive constant, so that (8.4) is the same as 
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(9.9) . If x,y e Q are such that \x-y\ < N^diy) then (9.9) follows directly from (8.13). Finally, if 
N~ l d(y) < \x-y\ < Nd(y), then, as before, d(x) < (1 +N)\x-y\ and 1 + is bounded from below 
and above by a positive constant. Then (8.5), (8.13), and (8.8) imply the bound on the left-hand 
side of (9.9) by C \x - y\ m ~i~\ a \ and hence, (9.9), using that d(y) « \x - y\. This finishes the proof 
of (9.4). 

Finally, we obtain from (9.4) the bound 

(9.10) \V m ^u(x)\ < C V fd(y) m -i- M 1-^-] \f a (y)\dy, x e Q. e > 0, 

Then, by the mapping properties of the Riesz potential we recover an estimate 

(9.11) ||V"-t M ||^ (£J) < C n Yj Hyy n ^ M+E M\ LP (a)' Prided that £ > 0, p> 



n-e 

\a\<m-j 



which, in turn, leads to (9.5). □ 
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